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COVARIANCE MATRIX OF THE POTTS MODEL 

Abstract. We consider the covariance matrix 

G"""(a; -y) = {q5{ax,m) q5{cjy,n)) - {qS{ax,m)) {qS{ay,n)) 

of the d-dimensional q'-states Potts model, rewriting it in terms of the connectivity, 
the finite-cluster connectivity and the infinite-cluster covariance in the random clus- 
ter representation of Fortuin and Kasteleyn. In any of the q ordered phases, we show 
that - in addition to the trivial eigenvalue - the matrix G"^'^(x — y) has one simple 
eigenvalue G\^-^^{x — y) and one {q — 2)-fold degenerate eigenvalue G^jj.(x — y). Fur- 
thermore, we identify the eigenvalues both in terms of representations of the unbroken 
symmetry group of the model, and in terms of connectivities and cluster covariances, 
thereby attributing algebraic significance to these stochastic geometric quantities. In 
addition to establishing the existence of the correlation lengths ^ . and ^ . corre- 
sponding to G^^lix - y) and G^^l{x - y), we show that ^^Vr(/3) > d'r(/3) ^r all 

(2) 

inverse temperatures f3. We also prove that ^^j^. coincides with the decay rate of the 
diameter of finite clusters. 

(2) 

For dimension d = 2 and q > 1, we establish a duality relation between ^ . and 
^freej the Correlation length of the two-point function with free boundary conditions: 
We show ^^jj.(/3) = 4^free(/3*) for all (3 > (3o, where (3* is the dual inverse temperature 
and f3o is the self-dual point. For systems with first-order transitions, this relation 
helps to resolve certain inconsistencies between recent exact and numerical work on 
correlation lengths at /3 = /3o- For systems with second order transitions, this relation 
implies Widom scaling. Namely, asssuming that ^free(/3) ~ \P — Pol^'^ as /3 t /3o, the 
duality relation gives i^l{(3*) ~ |/3* - /3o|"'^ as /3* j /3o with i> = u. 

In the course of proving the above results, we establish several properties of inde- 
pendent interest for the random cluster model, including left continuity of the inverse 

("2") 

correlation length l/Cfree(/3) and upper semicontinuity of the inverse length 1/^ . (/3) 
in all dimensions, and left continuity of the two-dimensional free boundary condition 
percolation probability at /3o- We also introduce DLR equations for the random 
cluster model and use them to establish ergodicity of the free measure. 

In order to prove the above results, we introduce a new class of events which we 
call decoupling events and two inequalities for these events. The first is similar to the 
FKG inequality, but holds for events which are neither increasing nor decreasing, and 
replaces independence in the standard percolation model; the second replaces the van 
den Berg - Kesten inequality. 



1. Introduction: Background and Discussion of Results 

The g-state Potts model has been the subject of increasing interest in recent 
years. On the one hand, it has been studied by probabihsts and statistical mechani- 
cists due to its relationship to the random cluster model (see [FK] and [ACCN]), 
where many of the known results for percolation are open and interesting problems. 
On the other hand, the phase transitions in the Potts model provide a paradigm 
for testing numerical methods developed for more complex transitions, such as de- 
confinement in lattice QCD: The Potts model is relatively easy to simulate with 
efficient algorithms (see e.g. [SW]), it can be tuned from a second-order through a 
weakly first-order to a strong first-order transition by varying the number of states 
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q, and many quantities of interest are explicitly known for dimension d = 2, thus 
allowing for a direct test of numerical methods. Finally, many of the exact re- 
sults on the Potts model have recently been shown to have fascinating algebraic 
interpretations (see e.g. section VII. B in the review of [W]). 

Motivated by discrepancies between recent exact and numerical results on the 
correlation length of the Potts model, we have undertaken to identify and study the 
relevant length scales in the problem. We relate these scales both to the algebraic 
structure of the unbroken symmetry group and to stochastic geometric quantities in 
the random cluster representation of the Potts model. In the process, we show that 
the some of the natural stochastic geometric quantities one defines in the random 
cluster representation - e.g. the finite-cluster connectivity - have independent al- 
gebraic significance. In two dimensions, we prove a relation between various scales 
which is an extension of known relations for percolation and the Ising magnet, and 
which establishes a strong form of Widom scaling for Potts models with continuous 
transitions. We also prove an analogue of this relation for two-dimensional Potts 
models with discontinuous transitions; this analogue helps to explain the apparent 
discrepancy between the exact and numerical results. 

Adopting a field theoretic perspective, we identify the relevant lengths in the 
model by studying the eigenvalues of the covariance matrix 

Grix -y) = {q5{a^, m) ; q5{ay, n)), . (1.1) 

Here ax G {0, ■■■ ,5 — 1} are the usual spins of the Potts model, 5(-, ■) is the 
Kronecker delta, ( ■ )^ is the expectation with respect to the infinite- volume state 
obtained from finite-volume states with "6" boundary conditions, and {A; B)^ = 
(AB)^ — {A)f^{B)ij is the truncated expectation of the functions A and B. 

In the disordered phase, we consider the covariance matrix with free boundary 
conditions, G^^^{x — y). We find that this is proportional to the standard two-point 
function, which in turn is equal to the connectivity function in the random cluster 
representation: 

G7^eeix-y) = {qS{m,n) -l){^-^{q6{ax,ay) -1))^^^^ = {q6{m,n) -l)Tiree{x-y) . 

(1.2) 
Here the connectivity, T{j-ee{x — y), is the probabililty with respect to the free bound- 
ary condition random cluster measure that x and y lie in the same component. That 
fhee{x — y) is equal to the two-point function in finite volume is well known both to 
probabilists and to numerical physicists, the latter of whom use this equivalence to 
measure the two-point function according to the "improved estimators" approach. 
Our only contribution here is to verify the equivalence in infinite volume. We note 
that, in the disordered phase, the covariance matrix contains no more information 
than the standard two-point function, or equivalently, the connectivity function. 

The problem is more subtle in the ordered phase, where we consider the matrix 
Q-mn^^ _ y^ with fixed constant boundary conditions, c G S = {0, ■ ■ ■ , g — 1}. 
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Defining the finite-cluster connectivity, T^fj-{x — y), to be the probability, in the so- 
called wired random cluster measure, that x and y lie in the same finite component, 
and the infinite- cluster covariance, Cwir(a; — y), to be the covariance, again in the 
wired measure, of the events that x and y lie in the infinite component, we prove 
that the matrix elements G'^^{x — y) are linear combinations of T^f^{x — y) and 
Cwiv{x-y), namely 

G^^{x -y) = {qd{m, n) - 1) r^U^ - y) + {qd{m, c) - l){qS{n, c) - 1)0^^^ - y) . 

(1.3) 

We remark that such a relation was not known previously, either in finite or infinite 
volume; furthermore, the proof of the infinite- volume limit involves some subtleties 
related to how the infinite cluster emerges from large finite clusters in the wired 
problem (for more details, see the remark following Proposition 3.4). 

Percolation analogues of T^^{x — y) and Cwir(a; — y) -in the absence of boundary 
conditions - have arisen previously in [CCGKS], where they appeared as a natu- 
ral decomposition of the truncated percolation connectivity in the ordered phase. 
There, however, they did not have independent signficance, appearing only as a 
sum. The question naturally arises whether they have independent significance 
here. Obviously, this is not the case for q = 2, for which (1.3) can be rewritten as 

Grix -y) = (25(m, n) - 1) {r^A^ - y) + C^.^x - y)) , 

involving again only the sum t^^x — y) + C^-^^x — y). 

For Q > 3, however, the fixed boundary condition covariance matrix G'^^{x — y) 
has a richer structure. We prove that it has a simple eigenvalue zero and a nontrivial 
simple eigenvalue 

G^2{x -y)=q T^Ax -y)+ q{q - l)C^Ux - y) , (1-4) 

both corresponding to the trival representation of the unbroken subgroup Sq-i of 
permutations of 5" \ {c}, as well as one {q — 2)-fold degenerate eigenvalue 

G^^l{x-y) = qT^Ax-y), (1.5) 

corresponding to the remaining orthogonal subspace.^ Thus we see that for Q > 3, 
the finite-cluster cluster connectivity, r^i^(x — y), has independent algebraic signifi- 
cance as an eigenvalue of the covariance matrix, and hence also physical significance 
in terms of the associated one-particle spectrum. As for the infinite cluster covari- 
ance Cwir(a^ ~ y)^ ^^ will show in Theorem 4.3 that its decay rate is equal to the 



^For the Ising model (ij = 2), G'^^{x — y) has only the trivial eigenvalue zero and the eigenvalue 

c'-'^hx-y). 

wir V ^z 
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decay rate of the eigenvalue G^j^. whenever the magnetization is positive. Thus al- 
though C^ir{x — y) does not have independent algebraic significance, its decay rate 
does. 

For completeness, we note that the free boundary condition matrix, Gf^^^{x — y), 
can be diagonalized as well, yielding a simple eigenvalue zero and a, {q — l)-fold 
degenerate eigenvalue 

Giree{x-y) = qXireeix -y) . (1.6) 

Given the eigenvalues (1.4)-(1.6), one naturally defines the inverse correlation 

lengths: 

1 1 

7 — 7m = ~ , ^™ VI logG'free(a;) , (1-7) 

Cfree(p) Ix^oo \x\ 

(1.8) 

and 

(1.9) 



1 


- lim ^logG'^lix 

\x\^oo X 

lim ^logGS(x) 

|x|— >CXD \X\ 


1 


£li(^) 



Here, as usual, (3 is the inverse temperature of the model. In all cases, the limits 
are taken so that x lies along a coordinate axis. In order to establish the existence 
of the limits, we return to the spin representation and use reflection positivity. 
We also give alternative subadditive proofs of the existence of the limits (1.7) and 
(1.9), which though more complicated than the reflection positivity arguments, have 
the advantage that they hold for non-integer q > 1 and can be used to establish 
additional properties. In particular, we use subadditivity to show left continuity 
of the inverse correlation length l/^free(/3) and upper semi continuity of the inverse 
correlation length l/^^in/^)- ^^ ^^^° ^^^ subadditivity arguments to prove that 
Wirl/^) is equal to several other geometrical correlation lengths in the problem, one 
of which is the decay rate of the diameter of finite clusters in the wired measure — 
a quantity which should be easiy accessible to numerical measurement. 

All three correlation lengths coincide in the high-temperature regime, where their 
common value is often denoted by ^dis(/3). In the low-temperature regime, we expect 
^free(/9) = oo. The uou-trivial correlation length in this regime is often denoted by 
^ord(/9)- Here, however, we see that for Q > 3, there are two a priori different non- 
trivial lengths, ^^jj.(/9) and $}^-^^{l3). Equations (1.4) and (1.5) immediately imply 
that 

&{^) > &i.^) , (1-10) 

so that the correlation length ^^j^. of the symmetric state (i.e. symmetric with re- 
spect to Sq-i) is not smaller than those of the unsymmetric states. An interesting 
open question is whether or not the inequality is strict. It is worth noting that in 
percolation, analogues of C^\^{x — y) and r^f^-ix — y) in the absence of boundary 
conditions have equal exponential decay rates [CCGKS], which here would imply 
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equality of ^^jj.(/3) and ^„ii.(/3)- However, it is not at all clear whether the Potts 
models for q > 3 should have analogous behavior. 

We return finally to our original question, namely the discrepancy between the 
exact and numerical correlation lengths of two-dimensional Potts models with dis- 
continuous transitions. Explicit calculations based on a mapping to the six-vertex 
model yielded a correlation length Cdisif^o) of the disordered phase at the self-dual 
point (3o ([BW], [KSZ], [Kl]) which disagreed with previous numerical measure- 
ments ([PL], [GI]) of the ordered correlation length at the transition point ^ord(/3o) 
by roughly a factor of 2, suggesting the possible relation ^ord(/9o) = |^dis(/3o) [BJ]. 
A continuous transition analogue of this relation is already known for both two- 
dimensional bond percolation, where ^{p) = |^(1 — p) has been rigorously es- 
tablished for all p > Pc [CCGKS], and the two-dimensional Ising magnet, where 
^{(3) = ^^{(3*) has been been established via exact solution for all (3 > (3o [MW]. 
Here, as usual, (3* is the dual inverse temperature. However, from our results dis- 
cussed above, we now know that the situation is more complicated in the g-state 
Potts model, Q > 3, than it is in either percolation or the Ising magnet, since in the 
ordered phase the Potts model has two a priori different correlation lengths. One 
of our principal results is a relation of the conjectured form in terms of the smaller 

(2) 

ordered correlation length, P .^. 

Our result follows from a dichotomy which we prove for all two-dimensional 
random cluster models with q > 1. In addition to the conjectured relation, the 
dichotomy implies Widom scaling for Potts models with continuous transitions. 
Let P^'^{(3) be the percolation probabilty in the free boundary condition random 
cluster measure. Our dichotomy is: If P^'^{(3*) = 0, then 

whereas if P^^^'^P*) > 0, then 

CireeW) = CHiP) = C^^iP) ■ (1-12) 

In order to interpret the dichotomy, we supplement it with the two-dimensional 
relation 

P--(/3)pfr-(r) = 0, (1.13) 

where P^^{f3) is the percolation probability in the wired measure, which is of 
course equal to the spontaneous magnetization M(/3). Note that (1.13) shows that 
P^^{(3*) > implies M{(3) = 0, so that (1.12) is simply the equality of the three 
correlation lengths in the high-temperature regime, as mentioned earlier. 

Our more interesting corollaries follow from the first branch of the dichotomy, i.e. 
the duality relation (1.11). In order to see this, we combine (1.13) with the obvious 
bound PZ'"i/3) > P^c^'^'^iP) to obtain pfr^^(/3) pfr<^^(/?*) = 0, so that P^'^(/?o) = 0. 
Since P^'^{(3*) is an increasing function of /9*, this in turn implies 

pfr^^(/3*) = forall /3>/?o. (1.14) 
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Equation (1.14) implies in particular that P^^{(3) is left continuous at the self- 
dual point Po- Moreover, it means that that the first branch of the dichotomy (i.e. 
equation (1.11)) holds throughout the low-temperature phase (3 > (3o- For systems 
with first-order transitions, this implies the conjectured relation at Po'- 

&(Po) = ICfreeiPo) . (1.15) 

For systems with second-order transitions, (1.11) is a generalization of the afore- 
mentioined results on two-dimensional percolation [CCGKS] and the Ising magnet 
[MW]. In particular, it gives a strong form of Widom scaling as /3 ^ (3o: If ^free(/3*) 
diverges with critical exponent z/, ^free(/3*) ~ |/3* — /3o|~^ as (3* t /3o, (l.H) implies 

(2) (2) 

that ^^jj.(/3) diverges with the same exponent: ^^jj.(/3) ~ \(3 — (3t\~'^ as /3 | /3o with 
u = u. 

As noted above, the interpretation (and in fact, the proof) of the dichotomy 
(1.11) and (1.12) requires the relation (1.13), which we obtain as a special case of 
a general two-dimensional result of Gandolfi, Keane and Russo [GKR]. However, in 
order to apply the [GKR] theorem, we need to know that the free random measure 
is ergodic, a result which we establish in all dimensions. We prove ergodicity by 
introducing suitable DLR equations ([D], [LR]) for the random cluster problem. 
Here the justification of the DLR equations is much more delicate than in standard 
spin systems due to the nonlocal nature of the random cluster weights: because of 
this nonlocality, the specification used to construct the DLR states is not quasilocal, 
and thus standard theorems do not apply. 

Before reviewing the organization of the paper, let us briefiy discuss our meth- 
ods. These methods are necessarily quite different from those used in the analysis of 
the Bernoulli percolation model, since the random cluster model lacks several prop- 
erties which are used extensively in percolation - namely, independence of events 
occurring on fixed disjoint sets and the van den Berg - Kesten [BK] inequality for 
events occurring on random disjoint sets. Moreover, the random cluster model has 
an additional feature - boundary conditions - which significantly complicates its 
analysis relative to the independent model. However, it is by actually focusing on 
the boundary conditions that we are able to circumvent the other difficulties and 
in fact derive two correlation inequalities which we expect will be useful in many 
other contexts. We do this by noting that in many cases, the events of interest 
carry with them boundary conditions which decouple them from other events and 
thus effectively overcome the coupling of the random cluster weights. This idea is 
formalized by introducing the notion of decoupling events. We use our decoupling 
events in formulating and proving two sets of inequalities which effectively replace 
independence and the BK inequality. The independence is replaced by a relation 
which resembles the FKG inequality, but contains two decoupling events and holds 
for a much larger class of events than the original FKG inequality - in particular, for 
events which are neither increasing nor decreasing. The BK inequality is replaced 
by a relation which resembles the independent BK inequality but contains a decou- 
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pling event. Both inequalities hold for any FKG measure, and thus in particular 
for the free and wired random cluster measures with q > 1. 

The organization of this paper is as follows. In first two parts of Section 2, we 
review the necessary properties of the standard spin and random cluster represen- 
tations of the Potts model. The third part of Section 2 contains our inequalities for 
decoupling events. In the last part of the section, we derive the DLR equation and 
establish ergodicity of the free measure. Section 3 is concerned with the covariance 
matrix. In the first part of the section, we derive the finite- and infinite- volume rep- 
resentations of the matrix with free and constant boundary conditions, in particular 
establishing the infinite- volume limits of rfree(a^ — y), ^wfri^ '~ v) ^'^d C^\-[{x — y) 
from their finite- volume analogues. The matrix is diagonalized in the second part of 

Section 3. Section 4 concerns the correlation lengths ^free, ^wir ^"^^ Wir- -'-^ ^^^ ^^^^ 
part of the section, we establish existence of the lengths using refiection positivity, 
as reviewed in the Appendix. The second part of the section concerns alternative 
characterizations of ^free, ^wir ^'^'^ Wir' pi'oved via subadditivity arguments and our 
inequalities for decoupling events. In particular, we show that l/^free is left continu- 
ous and l/^^ir is upper semicontinuous; we prove that ^^j^. is the decay rate of Cwir 
whenever the magnetization is positive; and we establish that ^^j^. coincides with 
the decay rate of the diameter of finite clusters, as well as with the limiting decay 
rate of connectivity functions for clusters in boxes. Section 5 contains our proof 
of the two-dimensional dichotomy (1.11) and (1.12), as well as derivations of a few 
results on two-dimensional percolation probabilities. The first part of this section 
contains a discussion of the heuristics of the duality relation (1.11) in terms of the 
behavior of interfaces in the system. In the second and third parts of the section, 
we prove upper and lower bounds of the form needed for the duality relation (1.11). 
Finally, in the fourth part of the section, we combine these upper and lower bounds 
with several results from Section 4 and the relation (1.13) to obtain our dichotomy. 

2. Preliminaries 

2.1. Definition of the Spin Model. 

We consider the g-states Potts ferromagnet, a model with spins a^ in the set 
5" = {0, 1, ■ ■ ■ , Q — 1}, Q > 1. In a finite volume A C Z^, the Hamiltonian with free 
boundary conditions is 

i^free(aA) = - Y. {K^^^^y)-^) (2.1) 

{x,y)eB{A) 

where the sum goes over the set -B(A) of all nearest neighbor pairs {x,y) for which 
both X and y lie in A. The Hamiltonian with c-boundary conditions, c & S = 
{0, 1,--- ,Q-1}, is 

H,{aA) = iffree((^A) " J^ (5((7„ Gy) - 1) , (2.2) 

aeaA 
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where dK = {x ^ A|dist(a:, A) = 1} is the (outer) boundary of A. Using the label b 
for "free" or c G 5", one introduces the partition function with boundary condition 
b as 

Z6(A) = ^e-^-f^''("'^) (2.3) 

where the sum runs over all configurations a\ : A ^ S,x ^-^ a^ and P is the inverse 
temperature (3 = l/ksT. 

As usual, an observable A with support suppA is a function A : cta ^ C which 
does not depend on the spin variables a^ for x ^ suppA. A local observable is an 
observable with a support suppA not depending on A. Expectation values of a local 
observable A are defined as 

(^)6,A = ^E^(^A)e-^^^(-). (2.4) 

If A and A are two local observables, one also considers the truncated expectation 
value defined by 

{A;A),^^ = {AA),^^-{A),^^{A),^^. (2.5) 

For observables of the form 

Ap = exp(z(cr,p)) = exp(z ^ a^p^) , 

X Esupp Ap 

where p is a function of finite support from Z^ into 5" = {0, 27r/g, ■ ■ • , 27r(g — !)/?}, 
the expectation values (^p)fi.ee a ^^^ monotone increasing (i.e. non-decreasing) in 
A, while the expectation values {Ap)^ ^ are monotone decreasing in A by Griffiths' 
second inequality [Grf], as generalized by Ginibre [Gi]. As a consequence, for b — 
"free" or 6 = 0, the thermodynamic limit 

(Ap), = lim (Ap),^^ (2.6) 

exists for all local observables of the form Ap = e'^^^'P' and hence for all local 
observables A. In (2.6), the limit may be taken through any increasing sequence of 
sets. Using the permutation symmetry of the Hamitonian (2.2), one concludes that 
the limit (2.6) exists for all boundary conditons b considered here (i.e. free or any 
constant boundary conditions). Also by Griffiths' second inequality ([Grf], [Gi]), 
the limit (2.6) is translation invariant. 

The order parameter of the Potts model is the magnetization 

M{P) = -^{q5{a,,,0)-l),, (2.7) 

where xq is an arbitrary point in Z'^ (recall that the infinite- volume states ( ■ )^ are 
translation invariant). It is known that M(/3) is increasing in /3 ([Gi], [ACCN]), 
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decreasing in q [ACCN], and that the infinite-volume states ( ■ )^, c G 5" are equal 
to ( ■ )fj.gg if and only if M(/3) = [ACCN]^. Defining the transition point 

Pt = inf{/3|M(/3) > 0}, (2.8) 

we remark that it is believed that M{Pt) is increasing in q, and that 

qc = maxJQ G N\M{(3t) = 0} (2.9) 

is 4 for (i = 2 and 2 for d > 2. The fact that M{(3t) > 0, i.e. the existence of a 
first-order phase transition, has been rigorously established for all d > 2 provided q 
is sufficiently large ([KoS], see also [LMR] and [LMMRS]). 

2.2. The Random Cluster Representation: Review of Basic Properties. 

It is often useful to reexpress the g-state Potts model as an integer value of a 
two-parameter interacting percolation model, the so-called random cluster model 
of Fortuin and Kasteleyn [FK]. In order to set our notation and state the results 
we will use in the rest of this paper, we briefiy review the derivation and some 
basic properties of the FK representation. The representation is defined in terms 
of configurations a; G O = {0,1}'^'', where B^ = {{x,y)\x,y G Z^} is the nearest- 
neighbor bond lattice. For subsets S C B^, the configuration space is denoted by 
Ob = {0,1}^. 

Let us start with the finite-volume partition funtion with free boundary condi- 
tions. We write the Gibbs factor e~^^f'^°°*^'^'^) as 



n 



e 



P{S{a^,ay) — l) 



{x,y}eB{A) 

and expand the product with the help of the identity 

e^^^^^'^'^'y'^-^^ = {l-p)+p6{a^,ay) where p = 1 - e''^ . (2.10) 

We identify each term of this expansion with a configuration u G Ob(a); uj is chosen 
so that it is zero on those bonds for which the factor in the product is p, and one 
on those bonds for which the factor is pd^axjO'y). Geometrically, we think of the 
bonds b = {x,y) for which uj{h) = 1 as occupied or ordered, and those for which 
u){h) — Q as vacant or disordered. With a slight abuse of notation, we sometimes 
use the symbol u to denote the set of occupied bonds in -B(A), and u'^ to denote 
the set of empty bonds in -B(A), see e.g. (2.11) below. 
Rewriting the Gibbs factor in expanded form, we obtain 

Zfree(A)= Y. H^^-vtV^^ n '^K,^^)- (2-11) 



■^Actually, [ACCN] proved that all infinite- volume Gibbs states are equal to ( • )fj.gg if and only 
if M(/3) = 0. 



10 C. BoRGS, J. T. Chayes 

Evaluating the sum over cta, we pick up a factor q for each connected component 
of the graph (A, cu) (regarding isolated points as separate clusters). Denoting the 
number of clusters in this graph by i^{uj), we find 

Zfree(A)= Yl (1-P)'"''P'"'9*^"^- (2.12) 

It is an easy exercise to generalize (2.12) to the expectations of local observables 
A = A{a). One obtains 

(^)free,A= E Giree,Ai^)EireeiA\u) (2.13) 

(^6nB(A) 

where 

Gfree,A(c^) = ,^-^il-py^V^^q*^^^ (2.14) 

Afreet, AJ 

is the weight of the configuration cu, while i?fi-ee("|'^) is an average over spins with 
the spins constrained to be constant on each connected cluster of to and with values 
for different clusters being chosen uniformly from {0, 1, ■ ■ • , g — 1}. We remark that 
for the purposes of interpreting expectations of this sort, it is often convenient to 
consider the joint distribution on the spin and bond variables with weights given by 
the terms in (2.11), as introduced implicitly in [SW] and explicitly in [ES]. In terms 
of this distribution, the expectation £'free('|'^) is an average over the conditional 
distribution of spins, given the bond variables. 

For constant boundary conditions, one obtains a similar representation, with the 
following differences (as noted in [ACCN]): 
i) The set -B(A) is replaced by the set i?+(A) of all nearest neighbor pairs {x, y) for 

which at least one of the two points x and y lies in A. 
ii) The points of the boundary dK are regarded as preconnected or wired, in the sense 

that these points are taken to be lying in one cluster. This of course modifies the 

value oi ^{uj). 
iii) The expectation E'free(^|'^) in (2.14) is replaced by Ec{A\uj), where the average 

is computed with the additional constraint that spins in clusters connected to 

the boundary now only assume the value a^ = c. 
We have: 

(^)c,A= Yl G^ir,Ai^)Ec{A\u;), (2.15) 

where 

Gwir,A(c^) = ^^Al-p)^^V^^q*^'^^ (2.16) 

and Zwir = Zlces ^c = qZo- 

We denote by //free,A(-) and //wir,A(-) the finite- volume measures defined by the 
weights (2.14) and (2.16), respectively. 
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Remark. The measures /Ufree.Al') and /Uwir,A(') are defined on the probabihty spaces 
(Ob(A), ^b(a)) and (fiB+(A), ^b+(a)), respectively. (In general, we use J^b to denote 
the cr-algebra generated by cylinder events A C O^.) It is sometimes convenient to 
extend these to measures on the full space (O, JF) by declaring all bonds in B>d\B{K) 
to be vacant for /Ufi-ee.Al')? and all bonds in B^ \ B^{K) to be occupied for /Uwii.,A(')- 
An important property of the FK representation is that it obeys the Harris- 
FKG inequality. This inequality, first proved for percolation in [H] and proved for 
a large class of models in [FKG], was established for the q > 1 random cluster 
representation in [F] (see also [ACCN]). We begin with the standard: 

Definition 2.1. Consider the natural partial order on bond configurations uj G O^, 
B C Brf, namely u ^ u' if a;(6) = 1 ^ ^'ip) = 1. A function / : 0^ ^ M is said to be 
increasing if it is nondecreasing with respect to this partial order, i.e. f{uj) < f{uj') 
for all LJ ~< Lj' . An event is said to be increasing if its indicator is an increasing 
function. Similarly, a function / is decreasing if the function — / is increasing, and 
an event is decreasing if its complement is increasing. 

A measure fx on (Ob, J^b) is said to be an FKG measure if it obeys the so-called 
Harris-FKG inequality 

lj{AinA2)>ij{Ai)f,{A2) (2.17) 

for all pairs of increasing events Ai, A2 G J-'b- It is said to be a strong FKG measure 
if for each cylinder event C G J-'b, the conditional measure n{- \ C) is an FKG 
measure. Finally, a measure fx on (Ob,JFb) is said to FKG dominate a measure u 
on (Ob, J-'b), denoted by z/ < (U, if i^{A) < n{A) for all increasing events A G J-'b- 

FKG 

Proposition 2.2. ([F], [ACCN]) Let q > I. Then the finite-volume free and wired 
FK measures, /Ufree,A cind /Uwir,A; o'^e strong FKG measures. 

Consequences (See e.g. [ACCN]). 
1) The finite- volume measures are monotonic in the volume: 

y"free,A < /Ufree,A' if A C A' (2.18) 

FKG 

and 

yUwir.A > /Uwir,A' if A C A' , (2.19) 

FKG 

from which it follows that the infinite-volume measures 

yUfree(-) = hm /Ufree,A(-) (2.20) 

A^Z^ 



and 



/Uwir(-) = lim /Uwir,A(-) (2-21) 



exist for all monotone local functions, and hence for all local functions. Further- 
more, these infinite-volume measures are translation invariant and inherit the 
strong FKG property. 
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2) The wired measures FKG dominate the free measures, i.e. 



and 



/^free,A < yUwir.A (2.22) 

FKG 



/^free S /-*wir • yZ.Zo) 

FKG 



Another useful property of these measures is that they have finite energy, a notion 
introduced by Newman and Schulman [NS]. 

Definition 2.3. Let B C B^, \B\ < oo, and (p G fls a configuration on B. If a; G O 
is a configuration on the fuU space, let (p{iv) be the configuration which agrees with 
(p on B and with uj on B'^: 

Finally, if A C O is an event, let (t){A) = {(p{u)) \ uj G A}. The measure /U on O is 
said to have finite energy if for every finite B cMd and for every (p G O^, 

^(A) >0 ^ /"(</>(^)) >0. 



It is easy to see that finite energy is equivalent to the statement: For each bond 
6, the conditional probability of the event that b is occupied, given the configuration 
on all the other bonds, is non-trivial: 

< fi{u{b) = l\u(b),b^b) < 1. 

For the free and wired measures, it was observed in [ACCN] that this probability 
can be explicitly calculated: 

{ if the endpoints of b are connected, 

regardless of ^(6) ^ ^2.24) 

p+gfi-p) Otherwise 

where fx = /Ufree or /Uwii-. Thus for all q' > 1 and all p 7^ 0, 1, the random cluster 
measures /Ufree and //wir have finite energy. Note that this is not true in all ran- 
dom cluster measures: Boundary conditions can impose constraints which exclude 
certain configurations. 

Given stationarity and finite energy, it follows immediately from a general result 
of Burton and Keane [BuK] that the infinite cluster is unique: 
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Proposition 2.4. For any q > 1 and any p G (0, 1), the free and wired random 
duster states have at most one infinite cluster with probability one. 

Since the Burton and Keane theorem requires only stationarity, it applies also to 
non-extremal states, and therefore allows the possibility of a convex combination of 
states with zero and one infinite cluster. If, in addition, the measures are ergodic, 
then at any given value of p, there is either zero or one infinite cluster with probabil- 
ity one. This is presumably the case for both the free and wired measures, although 
we only prove it for the free state (see subsection 2.4). Of course, ergodicity does 
not exclude the possibility that, for a fixed value of p, the wired state has an infinite 
cluster and the free state does not - indeed, for q large enough, this is exactly what 
happens at the transition point. 

2.3. Two Useful Inequalities. 

There are three main technical tools for factoring intersections of events in stan- 
dard Bernoulli percolation: the FKG inequality for monotone events, independence 
for events which occur on nonrandom disjoint sets, and the van den Berg - Kesten 
[BK] inequality for events which occur on random disjoint sets. As discussed in the 
last section, the free and wired random cluster measures obey an FKG inequality. 
However, due to the nonlocality of the weights (2.14) and (2.16), they satisfy nei- 
ther an independence condition nor a BK inequality. Indeed, it is clear from (2.24) 
that the probability of even a simple bond occupation event can be enhanced by 
the occurrence of some other event at an arbitrarily long distance from the bond in 
question. In this subsection, we provide alternatives to independence and the BK 
inequality for many events of interest in a general setting. 

As a substitute for independence of events occuring on nonrandom disjoint sets, 
we might try to use the FKG inequality as a bound, provided that the desired 
events are monotone. However, many of the events we care about - especially in 
the low-temperature phase - are not monotone. For example, the probability of a 
connection via finite clusters is the intersection of an increasing and a decreasing 
event. The presence of boundary conditions, which very often complicates proofs 
in the random cluster model, can be used to our advantage here. Certain boundary 
conditions decouple a set from its exterior. Many events of interest carry with them 
decoupling boundary conditions for the (random) sets on which they occur. We 
make this notion precise by introducing the definition of a decoupling event below. 
It turns out that, given this definition, it is possible to prove a general inequality 
which is similar to the FKG inequality and which replaces independence for events 
whose random boundaries occur within disjoint nonrandom sets. Our inequality 
holds for any FKG measure and for events which are intersections of arbitrary 
events with monotone decoupling events. 

As explained above, the BK inequality is certainly not true in general for the 
random cluster model - there are numerous examples in which the occurrence of 
one event enhances the occurrence of another. However, this enhancement cannot 
take place if the two events are decoupled from one another, in a sense to be made 
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precise in the definition below. Thus we prove a second inequahty, which replaces 
the BK inequality of Bernoulli percolation, and which holds for the intersection of 
an arbitrary event, an increasing event and a decreasing decoupling event. 

In Proposition 2.6 below, we actually present two versions of each of our inequal- 
ities: one which is easy to formulate (but not that useful), and a more involved one 
which is of the form needed for our applications. All of these inequalities hold for 
general FKG measures. We also give a useful corollary that concerns monotonicity 
in the volume and FKG domination in the random cluster model. We begin with 
the definition of a decoupling event. 

Definition 2.5. Given a probability space (O, JF, n) and events Ai, A2, -D G JF, we 
say that D is a decoupling event for Ai and A2, if 

fi{Ai nA2\D)= fi{Ai I D) fi{A2 I D) . (2.25) 

For brevity, we will sometimes say D decouples Ai from A2. 

While this definition makes sense in any probability space, it may be useful to 
illustrate it with a typical example from the random cluster model. Consider a set 
B cBd such that Md\B = BiU B2, Bi n B2 = 0. The event that the bonds of 
B are vacant then decouples any event Ai G J-'biUB from any event A2 G J-'b2UB- 
In this paper, such decoupling events typically occur when B is the boundary of a 
finite occupied cluster. Returning to the general context of Definition 2.5, we have: 

Proposition 2.6. Let {Q,J-',iJ,) be a probability space with O partially ordered and 
H an FKG measure with respect to this order. Then the following inequalities hold: 

1) The First Inequality 

i) Consider two arbitrary events Ai, A2 G J^ , and two increasing (or two decreas- 
ing) events Di, D2 E J-' such that Di decouples Ai from D2 while D2 decouples 
A2 from AiH Di. Then Ei = AiH Di and E'2 = ^2 H D2 obey the inequality 

fi{EinE2)>fi{Ei)fi{E2). (2.26) 

a) More generally, let Ei, i = 1, 2, be disjoint unions of the form 

E^= [j A,^k n A,fc, (2.27) 

keKi 

where Ki are countable index sets, Ai^k G JF are arbitrary events, Di^k G J-' are 
all increasing (or all decreasing) events, and Di^k decouples Ai^k from -D2,A;' while 
D2,k' decouples A2^k' from Ai^k n-Di,A: for all k G Ki and k' G K2. Then Ei and 
E2 obey the inequality (2.26). 

2) The Second Inequality 

i) Let Ai E J-' be an increasing event, A2 E J-' be arbitrary, and D E J-' be a 
decreasing event which decouples Ai from A2. Then 

^i{Ai n L> n A2) < ii{Ai)ii{D n A2) < ii{Ai)ii{A2). (2.28) 
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ii) More generally, let Ai E J-' be an increasing event, and let A2 E J-' and D E J-' 
be events for which DnA2 can be rewritten as a disjoint union of the form (2.27), 
with D2,k decreasing events that decouple Ai from A2,k for all k G K2. Then the 
bound (2.28) remains valid. 

Proof. Rewriting the left hand side of (2.26) as 

^{D2) fi{Ai n Di n A2 I D2) 
and using the fact that D2 decouples A2 from AiH Di, we obtain 

//(Ai n Di n A2 n D2) = fi{Ai n Di n D2) fi{A2 \ D2). 

Applying the same procedure to the term |u(Ai nDi \~\D2) and using the decoupling 
event Di, we get 

At(Ai n Di n A2 n D2) = ij{Di n D2) n{Ai \ Di) ii{A2 \ D2) 

which by the FKG inequality (2.17) implies (2.26). Part lii) of the proposition then 
follows from the countable additivity of the measure ^ and the fact that the events 
El and E2 are disjoint unions of events for which (2.26) is valid. 
In order to prove 2i), we observe that 

Ai(Ai n L> n A2) = ^i{D) ^i{Ai I D) ^i{A2 I D) 

by the definition of conditional expectations and (2.25). Using the FKG inequality 
(2.17) to bound ^{Ai \ D) by //(Ai), the bound (2.28) now follows. Again, 2ii) 
follows from 2i) and the countable additivity of the measure. D 

Remark. It is clear from the above proof that the inequality (2.26) is reversed if one 
of the two decoupling events Di and D2 is increasing and the other is decreasing. 
Similarly, the first inequality in (2.28) is reversed if Ai and D are both decreasing 
or both increasing. 

Corollary. Let q > 1, h. d U^ and b = "wir" or "free". Consider the random 
cluster measure Hb,A and the corresponding probability space (O^j^, J-'Bt^ A*&,a); where 
Bb = B+{A) ifb= "wzr" and B^ = S(A) if b = "free". Let B C Bf,, and let E be 
an event of the form (2.27) where the index set is the set of all subsets of B, i.e. 

E= (j AsDDs, 

SCB 

with As G J-'bi, arbitrary events, and Ds G J-'sb decreasing events that decouple As 
from all events in J-'b,,\b- If the events Ds are decoupling events with respect to the 
measure //free, A; then 

l^iree,A'iE)> IJ'ireeM^) providcd A' C A and B C B{A'). (2.29) 
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If the events Ds are decoupling events with respect to the measure /Uwir,A; then 

IJ'wir,A'{E) < iJ^ir,A{E) provided A' C A and S C S+(A') (2.30) 

and 

l^whME) < l^heeM^) provided B C B{A). (2.31) 

Proof. Let Ai be the event that aU bonds in -B(A) \ -B(A') are vacant. Then Ai G 
^B{A)\B{A') C J-'b{A)\b is decoupled from As by the event Ds- Since both Ai and 
Ds are decreasing events, 

Hbee,A{^ I ^l) > /Ufree,A(-£^) 

by the remark following the proof of Proposition 2.6. Observing that Hhee,A'{E) = 
fJ'free,A{E \ Ai), this proves (2.29). Defining Ai as the event that all bonds in 
i?+(A) \i?+(A') are occupied (all bonds in i?(A)+ \i?(A) are vacant), the remaining 
two inequalities are proved in the same way. D 

In order to illustrate the utility of Proposition 2.6, we conclude this subsection 
with applications of each of the two inequalities. These applications will be needed 
in our subsequent analysis and may be of independent interest. As usual, we denote 
by C{x) = C{x; uj) the set of occupied bonds connected to x in the configuration a;, 
and define {x ^^ y} as the event that x is connected to y by a finite path of occupied 
bonds. We also define, for each finite set A C Z'^ and any two points x, y G A, the 
event i?^^(A) that x and y are connected by a cluster C{x) C B{A). 

Proposition 2.7. Let q > 1, A C Z'^ be finite or infinite, and let fx = /Uwir.A or 
/Ufree,A. Then for all finite Ai, A2 C A with B+{Ai)nB{A2) = S(Ai) nS+(A2)'= 0, 
and all x,y E Ai, z,w E A2, 

l^iRt^^yiA,) n R%iA2)) > fiiRl^iA^)) ix{R%{A2)). 

Proposition 2.8. Let q > 1, A d U^ he finite, and let n = /Uwir,A or /ifree,A- Let 
X, y, z,w E A, and let D be the event {x ^ z} H {y ^ w}. Then 

n{{x ^^ y} n D n {z ^^ w}) < n{x ^^ y) n{z ^^ w). 

Clearly, Proposition 2.7 is an application of the first inequality in Proposition 
2.6 - the connections in question occur on fixed disjoint sets, -B(Ai) and i?(A2), 
and due to the finiteness of the clusters, each connection carries its own decoupling 
event. Note that if B~^{Ai) fl i?"'"(A2) = 0, then in percolation, the probability of 
the intersection of the events in Proposition 2.7 would factor exactly. Here our 
first inequality replaces this independence. In fact, given that the decoupling events 
can overlap. Proposition 2.7 gives a new result even in the case of percolation. 
Proposition 2.8 is an application of the second inequality in Proposition 2.6 - the 
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connections in question occur on random disjoint sets separated by the decoupling 
event D. This obviously replaces the BK inequality. Note that, in marked contrast 
to percolation, the inequality would fail to hold if we removed the decoupling event. 

Proof of Proposition 2.7. Introducing Bi as the family of all sets B C -B(Ai) such 
that B connects x to y, and B2 as the family of all B C -B(A2) connecting z to w, 
we decompose i?^^(Ai) and i?f^(A2) as 

Rl^y{K,)= \J{C{x)=B}= \J{oOB = l}n{u9,B = Q} 
BeBi BeBi 



and 



R':%.i^2) = U {C{Z) =B}= [j{uB = l}n {U9*B = 0} 
BeHa BeBi 



Here ujb is the configuration uj restricted to the set B and d^B is the set of all 
bonds in Bd \ S which are connected to B. Observing that for all B & Bi and 
all B e B2, DiB = {uJd*B = 0} decouples Ai^b = Wb = 1} from all events in 
JF^c D ^B+(A2), while D2 B = {^Q*B = 0} decouples A^ j^ = {cu^ = 1} from all 
events in JF^^ d !Fb+{Kx)^ ^'^^ easily verifies that i?^^(Ai) and i?f^(A2) are events 
of the form considered in part lii) of Proposition 2.6. D 

Proof of Proposition 2.8. Defining Ai = {x ^^ y} and A2 = {z ^^ w}, we rewrite 
Ai as the disjoint union Af " U A^"^, with 

Af" = Ain{x^ dA} 

and 

A';"^ = Ain{x^ dA} . 

Notice that //free,A(^i"0 = Oj since with free boundary conditions, x cannot be 
connected to the outer boundary dA = {x ^ A|dist(x, A) = 1}. Introducing the 
family Bi of sets B C B{A) that connect x to y but do not connect x to 2; or y to 
w, we then decompose Ai H D a,s 

Af n I) = U {uB = 1} n {ujd*B = 0} . 

BeBi 

Observing that for all B E Bi, the event {uJd*B = 0} decouples A2 from the event 
{ujb = 1}, we obtain 

/.(Af nDnA2)< //(Af- n D) ^(^2) < ^(Af ) ^(^2) , 

where we have used the second inequality of Proposition 2.6 in the first step. This 
completes the proof for the free measure. 
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In order to complete the proof for the wired measure, we wiU show 

^{Af nDnA2)< i^{Af) ij{A2) . 

To this end, we define 

Af" = A2n{z^ dA} . 

Since the wiring would connect two points if they were both connected to the bound- 
ary, we have 

A'f n D n A2 = A'f n L> n Af '^ 

with probability one with respect to the wired measure. Applying the same strategy 
as before, we then obtain 

^{^f r\Dr\A2)= ^i{^f non A^"") < fi{Af) //(Af) < ^{Af) //(A2) , 

as claimed. D 

Remarks. 

1) As can be seen from the above proof, the finite- volume free measure actually 
obeys the stronger inequality 

/^free,A({a^ ^ V} H D H {z ^ w}) < ^free,A{{^ ^ V} ^ D) /Ufree,A(^ ^ w) . 

2) Using uniqueness of the infinite cluster [BuK] (see also Proposition 2.4 above), 
it follows immediately that A^"^ n D n A2 = A^"^ nDn Af " with probability one 
with respect to the infinite- volume measures //free and /Uwir- Hence Proposition 
2.8 holds for these measures as well. 

2.4. DLR Equations and States of the Random Cluster Model. 

In this subsection, we introduce the notion of (unconstrained) DLR states for the 
random cluster model, prove that the free measure is such a state, and use this to 
show that it is ergodic — a property we will need in our subsequent analysis. It is 
usually straightforward to establish such results by invoking the general theory of 
Gibbs states (see e.g. [Pr] and [Ge]). However, the general theory requires that the 
finite- volume expectations used to construct the DLR states are quasilocal functions 
of the boundary conditions, a property which fails to hold here due to the nonlo- 
cality of the random cluster weights. Thus the DLR equation has to be established 
explicitly. 

We start by defining finite- volume measures with general unconstrained boundary 
conditions — conditions which permit any component to be connected to any other 
component. The set of states generated by all such boundary conditions is quite 
natural in the random cluster model. A larger class including constrained states 
will be discussed briefiy at the end of this subsection. Each measure is defined on 
an arbitrary finite set of bonds B dMa with boundary 

dB = {x eU^\ 3y, zeU^ with (x, y) e B , {x, z) G B^} . 
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We specify the boundary condition by introducing a wiring diagram, W, which is a 
disjoint partition of dB into nw = 1, ■ ■ ■ , \dB\ components: 

nw 

w = {Wi,--- ,Wn^} with dB=[jw,, w,nWj=(l} ifiy^j. 

i=l 

We denote by yV{dB) the set of all such wiring diagrams - i.e. the set of all disjoint 
partitions of dB. Each component, Wi, of the wiring diagram W is considered to 
be preconnected or wired, so that all bonds b E B connected to points of Wi are 
regarded as being connected to each other. The number of components ^(cu) is 
then computed as usual. The random cluster weight 

GwM^) = ^^A'^-p)^''V'^^q*^^^ (2.32) 

defines the finite-volume measure hw,b{-)- Denoting by Wf^ee the partition with 
nw = \dB\ components and by VFwir the partition with only a single component, 
we see that HireeM') = I^Wt,,^,,B{A)i-) and /Uwir,A(-) = IJ'W„i,,B+{A)i-)i so that the 
free and (fully) wired measures are just special cases of nw,Bi-)- Note that among 
the measures ^w,b{') are some that cannot be obtained as transforms of any finite- 
volume states in the spin system, namely those in which W has more than q com- 
ponents Wi with |I4^i| > 2. 

There is a natural partial order on the set W{dB). If W, W G W{dB), we say 
that W is coarser than W , denoted by W >~ W , \i for each W^ G W there exist 
Wi^,Wi^,--- , Wi^ G W such that W[ = Uf^iWi.. Notice that Wf^ee is the least 
coarse and VFwir is the most coarse of all wiring diagrams. Moreover if W' >- W, 
then Hw',B dominates Hw,b in the sense of FKG (see Definition 2.1 above). 

Each configuration w G O induces a wiring diagram on each finite set B G M^. 
The induced wiring diagram W{B, uj) is a partition into components of dB, each of 
which is connected using occupied bonds in cusc. Thus each a; G O gives rise to a 
sequence of induced finite- volume measures Hw{b,u}),b for any increasing sequence 
of sets B C Brf. Henceforth we will extend the induced finite-volume measure 
Hw{B,u}),B to a measure on the full space (0,JF) by declaring all bonds in B'^ to 
have the configuration specified by uj. (Compare this to our extensions of of HfTee,A 
and /Uwir,A discussed in the remark following equation (2.16).) Using the form 
(2.32) of the weights Gw,b and our definition of induced wiring diagrams, it is 
straightforward to check that the (extended) induced finite-volume measures obey 
the consistency condition 

IJ'WiB,uj),B{A) = / IJ'WiB,Lu),B{duj) ^J'W{B,a,),Bi'^) (2.33) 

for all local events A G JF, any finite set B and all S C -B. 

For each finite B, we may define the function ttb '■ (^, ^) ^ M by nB{A\u) = 
I^w{b,uj),b{^)- Since the family 7 = {ttb I B C B^, |-B| < 00} is a set of proper 



20 C. BoRGS, J. T. Chayes 

probability kernels obeying the consistency condition (2.33), 7 is a specification in 
the sense of [Pr]. 

A DLR equation ([D], [LR]) is just an infinite- volume analogue of a consistency 
condition like (2.33). Thus we introduce the (unconstrained) DLR equation for an 
infinite-volume random cluster state fu.: 

fiiA) = / ^{du) ^w{b,co),b{A) (2.34) 

where A G JF is any local observable and B cMdis any finite set. As usual, the DLR 
equation (2.34) — if it holds — allows us to write the infinite-volume expectation 
of A as an average over finite- volume expectations. It is closed in the sense that the 
average is computed with respect to the given measure //. Note that this is different 
from the equation for states given in [ACCN], where a random cluster measure 
was obtained as a transform of a measure obeying the DLR equation in the spin 
system. On the other hand, a DLR equation was implicit in the discussion of states 
in [Grm]; there, however, the question of existence of solutions to the equation was 
not addressed. 

Let us denote the set of states obeying (2.34) by ^ = ^(7), where as above 7 
denotes the specification. States n E Q will be called DLR states or Gibbs states. 
A priori it is not clear whether Q is nonempty, i.e. whether there exists any n 
satsifying (2.34). One might try to construct such a /U as a subsequential limit 
of finite-volume measures hw,b — which clearly exists by compactness — but the 
question of whether such a limit obeys (2.34) involves a delicate interchange of 
limits. The theory of Gibbs states ([Pr], [Ge]) provides general conditions under 
which (2.34) is satisfied, one of which is quasilocality of the specification. 

A function / is quasilocal if it can be approximated in the supremum norm by 
local functions, a property which is equivalent ([Ge], Remark 2.21) to the statement 

sup I f{u) - f{Ti) HO as S ^ Bd . 

uj,ri: ujB=riB 

A specification {ttb} is quasilocal if the functions T^BiA,-) are quasilocal for all 
finite S C Bd and all local events A G JF. 

Unfortunately, due to nonlocality of the weights Gw,b, our specification is not 
quasilocal. For example, the probability of the simple event {uj{h) = 1}, conditioned 
on the bonds in B^^ \ {6}, changes discontinuously depending on whether or not the 
endpoints of b are connected by a path (of any length) in B^ \ {b} (see equation 
(2.24)). The general theory of Gibbs states can therefore not be applied here. 
However, we can verify the DLR equation (2.34) explicitly in the case of the free 
measure: 

Proposition 2.9. For all q > 1 and < (3 < 00, /Ufree ^ Q- 

Proof. Let S C B^ be a finite set and A E J-' a local event. We wish to show 

l^heeiA) = / IJheeiduj) IJ^W{B,u;),b{A) . (2.35) 



COVARIANCE MATRIX OF THE POTTS MODEL 21 

By the finite-volume consistency condition (2.33) and convergence of the finite- 
volume measures (2.20), it suflaces to prove 

liin jJfreeMd^) IJ'W{B,u;),b{A) = IJ^heeidiv) IJ'W{B,u;),b{A) . (2.36) 

^^^dj J 

Inserting the partition of unity J2weW(dB) ^{W{b,uj)=w} = 1 iiito (2.36) and noting 
that hw,b{A) is independent of A, we see that it is enough to prove 

lim ^i,ee,Ai{W{B, u) = W}) = l^i,ee{{W{B, u) = W}) , (2.37) 

i.e. that the probability of a given wiring diagram converges. 

Let RwiiB'^) denote the event that all sites within the set Wi are connected to 
each other via bonds in S'^, let SwiB'^) = Hvk eVK -^W'i(-^'^)' ^^'^ ^^^ Nwi,Wj{B'^) 
denote the event that none of the sites in Wi is connected to any of the sites in Wj 
via bonds in B^. Then 

{W{B,u;) = W} = Sw{B^) n f| Nw,,w,{B') ■ (2.38) 

By inclusion-exclusion, it is not hard to show that if -B(A) D B, then 

^l,reeA{W{B,Uj) = W}) = J^ k^{W) ^free,A(%(S^)) , (2.39) 

where the sum is over W coarser than W (see the definition a paragraph below 
(2.32)), and ky^{W) G Z are computable coefficients with kw{W) = 1. Thus by 
(2.37)-(2.39), we only need to show that for aU W G W{dB) 

hni ^free,A(%(i?^)) = ^free(%(i?^)) • (2.40) 

A ^i^d 

Let W G W{dB) and choose A such that -B(A) D B. Due to the free boundary 
conditions on dA, the argument of the left hand side of (2.40) can be rewritten as 

^free,A(^VK(5')) = ^free,A(^VK (5(^) \ B)). (2.41) 

Approximating the wiring event S^r{B'^) by local events, we see that the right hand 
side of (2.40) is actually a double limit: 

= hm lim ^free,A(%(S(A')\S)). (2.42) 

A' ^&d A -^&d 
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Thus by (2.40)-(2.42), we must show 

hm ^free,A(^w.(S(A) \ B)) = Hm hm fn^^^^SwiBi^') \ B)) • (2.43) 

A -^£d A' ^/,d A ^&d 

In order to prove this, we note that for all A' C A 

l^iree,A'iSwiBiA')\B)) < ^free,A(%(S(A') \S)) < ^free,A(^w^(S(A) \ S)) , (2.44) 

where the first inequality follows from the monotonicity property (2.18) and the 
second is just a consequence of Syy{B{A') \B) C S^^{B(A) \ B) if A' C A. Taking 
the hmits A ^ Z^ and A' ^ Za, equation (2.44) yields (2.43) and hence (2.35). D 

Remarks. 

1) The only property of the free measure that was used to reduce the proposition 
to equation (2.40) was convergence of the finite-volume measures. Thus the 
wired analogue of (2.40) — i.e. convergence of the probability of the wiring 
events Sy^{B'^) with respect to the finite- volume wired measures — is sufficient 
to prove //wir £ Q- Unfortunately, however, //„ir does not have nice monotonicity 
properties like those in equation (2.44). 

2) Equation (2.43) of the proof is our first example of the problem of interchange 
of limits which arises again in Proposition 3.4 and in many theorems in Section 
4. Whenever we deal with the infinite-volume limit of an event which is not 
confined to a finite volume, we encounter a double limit — one for the construc- 
tion of the infinite-volume measure and the other for the approximation of the 
given event by local events. Hence the problem of interchange of limits. This 
problem does not arise in percolation because the measure is defined directly in 
the infinite-volume limit. Here, when we can deal with interchange, it is usually 
accomplished via either simple FKG monotonicity (equations (2.18) and (2.19)) 
or our monotonicity involving decoupling events (corollary to Proposition 2.6). 

It is now straightforward to show that /Ufree is ergodic. We have: 

Theorem 2.10. Let H he any nontrivial subgroup of the translation group and 
let Qo '^ G be the set of all H -invariant DLR states. Then for all q > 1, ii{,:ee is 
extremal in Qo and hence is H-ergodic. 

Proof. As noted earlier, VFfree is the least coarse of all wiring diagrams, so that 

/^w^f_,B < m,B for all W e W{dB) , (2.45) 

FKG 

and thus by convergence of the measure (2.20) 

A^free < A* for all |U G ^ . (2.46) 

FKG 

Given that //free £ G (Proposition 2.9), it follows immediately from (2.46) that 
/Ufree is extremal in Q and hence also in Qo (since /Ufree is of course i7-invariant). 
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Ergodicity then follows from the fact that all extremal measures in Q^ are -ff-ergodic 
([Pr], Theorem 4.1). D 

Remarks. 

1) FKG Ordering of States: Using the fact that the wired state is the coarsest of all 
states, we have analogues of (2.45) and (2.46) for the wired measure, and thus 

A^free < A* < A^wir for all fi E G ■ (2.47) 

FKG FKG 

Note of course that this does not imply n^ij. G Q. 

2) The Size ofQ: By Proposition 2.9, /Ufree G ^ so that \Q\ > 1 for all g > 1 and all in- 
verse temperatures (3. Let P^'^i^f]) denote the percolation probabilty in the wired 
measure, which of course coincides with the magnetization for integer q. Accord- 
ing to a result of [ACCN] (Theorem A.2), whenever P^^'(/3) = (i.e. (3 < Pt for 
systems with second-order transitions and jS < fit for those with first-order tran- 
sitions) /Ufree = A^wir, SO that by (2.47) and Proposition 2.9, \Q\ — 1. It is expected 
that 1^1 = 1 also for (3 > f3t, but there are only incomplete results ior d = 2: The 
two-dimensional dual of the [ACCN] result says P^^^{/3*) = implies \Q\ = 1, i.e. 
there is one state for (3 > (3^, which presumably coincides with (3t (see also [Grm]). 
However, one expects more states at the transition point in systems with first- 
order transitions. For q large enough and d = 2, convergent expansions ([KoS], 
[LMR]) can be used to show that there are q + 1 distinct translation-invariant 
spin states (which transform into two distinct translation-invariant random clus- 
ter states - the free and the wired). There are presumably no non-translation- 
invariant states. Thus we expect \Q\ = 2 for (3 = (3t and q large enough in d = 2. 
In three dimensions, convergent expansions [MMRS] can be used to show that for 
q large enough, in addition to the translation-invariant states discussed above, 
there are infinitely many non-translation-invariant "Dobrushin-type" states cor- 
responding here to states constructed from wiring diagrams which coincide with 
Wwir above a certain hyperplane and with Wfree below that plane. We expect that 
these expansions can also be used to show that these non-translation-invariant 
states satisfy our DLR equation (2.34), so that aX (3 = (3t, \Q\ = oo for g large 
enough in (i > 3, in contrast to the conjecture of [Grm]. 

3) States with Constraints: In the remark above, we mentioned "Dobrushin-type" 
states which we expect to be in ^; these states were constructed from a combina- 
tion of wired and free boundary conditions. There are, however, many Dobrushin- 
type states in the spin system whose transforms are not in ^ — namely, mixed 
states in which various components of the boundary have different values of the 
spin. In the random cluster model, these correspond to states with constraints 
— certain components cannot be connected to other components. Therefore, in 
order to formulate DLR equations for these states, one has to supplement our 
wiring diagrams with some notion of constraints. While this is possible for in- 
dividual finite-volume states, it is not clear how constraints should be induced 
by a given configuration a; G O, nor whether the resulting measures would obey 
even finite- volume consistency conditions. 
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3. The Covariance Matrix 

3.1. The Random Cluster Representation of the Covariance Matrix. 

In this section, we rewrite the covariance matrices with free and constant bound- 
ary conditions, 

Gt:,{x -y) = {q5ia,, m) ; q5{ay, n))^,,, (3.1) 

and 

G'^^ix -y) = {qd{a,, m) ; qd{ay, n))^ , (3.2) 

in terms of the random cluster representation of Fortuin and Kasteleyn [FK]. We 
do this by first deriving finite-volume expressions and then taking infinite-volume 
limits. 

3.1.1. The Covariance Matrix in Finite Volume. 

Before deriving our representation for the covariance matrix, we recall the corre- 
sponding result for the (finite-volume) magnetization 

M,(/3,A) = J—{qS{a,,0)-l)^j,. (3.3) 

q — 1 

Using the symbol X ^^ Y for the event that the set X is connected to the set Y by 
a finite path of occupied bonds, the expression (2.15) almost immediately gives 

M,(/?, A) = //wir,A(a: ^ dA) . (3.4) 

For future reference, we note that this can be easily generalized to the expectation 
of e'P"- with pe S\{0} = {27v/q, ■ ■ ■2Ti{q - l)/q}. We obtain 

(e'''"^)o,A = /«wir,A(x^c>A) if peS\{Q}. (3.5) 

We begin by considering the finite- volume two-point function with free boundary 
conditions, 

G'SeeA^^ y) = {(iK^x, to) ; q5{ay, n))^^^^^^ . (3.6) 

Using the fact that {q5{ax,m))^^^^ j^ = 1 for all m and all a; G A, we first rewrite 
^free a(^' y) ^s ^^ untruncated expectation value 

G^:e,A(^,y) = {{q6{a^,m) - l){qS{ay,n) - 1))^^,,^^. (3.7) 

Now observe that £'free((<?5(o"x7 to) — l){qd{ayjn) — l)\uj) = if x and y are not 
connected in the configuration cu, while -Efree((?5(a"x, to) — l){qd{ay,n) — l)\uj) = 
q5{m, n) — 1 if X and y are connected. Thus, defining the connectivity in the FK 
representation 

Tfree,A(a;, ?/) = /Ufree,A(a^ ^ ?/) , (3.8) 

we obtain the following: 
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Lemma 3.1. The finite-volume covariance matrix with free boundary conditions 
has the representation 

<^w,a(3^. y) = iQ^{m, n) - l)rfree,A(a:^, y) ■ (3.9) 

Remark. The result (2.9) in [ACCN] for the usual two-point function, 

-—^{qSia^^ay) - l)free,A = ^reeM^^y) ' 

is proportional to the trace of our expression (3.9). 

Next, we rewrite the finite- volume covariance matrix with constant boundary 
conditions, 

GT,a(^ -y) = (?5(a,,m) ; q5{ay,n))^^^. (3.10) 

To this end, we define the finite-cluster connectivity 

r^,^^{x, y) = ;Uwir,A({x ^y}n{x^ M}), (3.11) 

and the covariance of the events that x and y are connected to the boundary dA. 

Cwir,A(a^, y) = /Uwir,A({a^ ^ dA} n{y ^ dA)} - /Uwir,A(a^ ^ M);Uwir,A(y ^ dA). 

(3.12) 

We have: 

Lemma 3.2. The finite-volume covariance matrix with constant boundary condi- 
tions has the representation 

GT,l{x, y) = {qd{m, n) - 1) t^,^^{x, y) + {qd{m, c) - l){q5{n, c) - l)C„ir,A(x, y) . 

(3.13) 



Proof. By the symmetry of the model, it is enough to establish the lemma for c = 0. 
In a first step, we prove a similar relation for (e*^*^^ ; e~*^ '^^)q ^, namely 



/g^P<T. 



; e-'P'''y)o,A = (1 - 5(p, 0))(1 - 5{p\ 0)) (Cwir,A(x, y) + 5(p, p')tXk{x. v)) ■ 

(3.14) 



Assume w.l.o.g. that p 7^ and p' 7^ 0, since otherwise (e*^'^== ; e *^'^^)oa ~ 0. 
Then recalling the definition of truncated expectation values and observing that by 

(3.5) 

(e'''"^)o,A(e"'^'"^)o,A = /«wir,A(2: ^ dA)^i^,,^p,{y ^ M), 
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the proof of (3.14) reduces to showing that 

^gW.g-^pV,^^^^ = ^wir,A({x ^ dA} n{y^ dA}) + S{p,p')rXA{^, v) ■ (3-15) 

We consider the cases p — p' and p ^ p' separately: If p 7^ p', the expectation 
i?o(e*^'^"=e~*^ '^y\uj) is zero unless both x and y are connected to the boundary, in 
which case £'o(e*^'^'"e~'^ "y\uj) = 1 . As a consequence, 

(e^^"^e-^p''^-)o,, = ;Uwir,A({x^aA}n{y^aA}) if p ^ p' . (3.16) 

If p = p', we consider two cases: either x <-> dA in the configuration a; or x jA 9A. In 
the first case, Eo{e'P''-e-^P'''y\uj) = 1 if y ^ ^A as weU, and Eo{e'P''-e-^P'''y\uj) = 
if y jA 9A, yielding a contribution of /Uwir,A(a^ ^^ (?A and y ^^ dA). In the second 
case, £'o(e'P'^^e~*^ ""''la;) = 1 if x ^ y, and £'o(e'P'^^e~*^ ""''la;) = if x jA y, yielding 
r^Ai^^y)- Thus 

(e^^"^e-'P'<^-)o^^ = ^„ir,A({x--aA}n{y^aA})+rfirr,A(x,2/) if p = p'. (3.17) 

Equations (3.16) and (3.17) establish (3.15) and hence (3.14). 

Given (3.14), the proof of the lemma is an easy exercise: observing that the 
delta functions q5{aj;,m) and q6{ay,n) can be rewritten as ^pg^ e'^*^'^'"""^^ and 
X]p'65^~*^ (<^y-"-)^ respectively, we multiply both sides of (3.14) by e'*^^ n-pm) ^^^ 
sum over p and p' to obtain (3.13) for c = 0. D 

3.1.2. The Covariance Matrix in Infinite Volume. 

In this subsection, we extend our representations of the covariance matrix with 
free and wired boundary conditions to the infinite volume. To this end, we again 
denote by C{x) = C{x;uj) the set of occupied bonds connected to x in the configu- 
ration u), and define the (translation-invariant) analogues of expression (3.8) for the 
connectivity, 

Theeix -y) ^ IJ^freeix ^ y) , (3.18) 

expression (3.11) for the finite- cluster connectivity, 

r^"j.(x -y) = /Uwir(a; ^ y and \C{x)\ < oo) , (3.19) 

and expression (3.12) for the covariance, 

Cwir(x -y) = Cov^_(|C(x)| = cx), \C{y)\ = oo) . (3.20) 

where in general Cov^(A, B) = n{A n B) — n{A)n{B) is the covariance of events A 
and B with respect to a measure //. We call the function C^hix—y) defined in (3.20) 
the infinite- cluster covariance. Our infinite-volume representation is contained in: 
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Theorem 3.3. The covariance matrices G'^^^{x — y) and G'^^{x — y) can be ex- 
pressed as 

GTe(^-y) = (qSi^.n) - l)r(x,y) (3.21) 

and 

G^^{x -y) = {q5{m, n) - 1) t^^{x - y) + {q5{m, c) - l){q5{n, c) - l)C^-Ax - y) • 

(3.22) 



Proof. Given the corresponding finite- volume statements in Lemma 3.1 and Lemma 
3.2, the theorem is an immediate consequence of the following proposition. 

Proposition 3.4. Let q > 1 be real, let Tiree,A{x,y), T^f^j^{x,y) and Cwir,A(ic, y) 
be the quantities defined in equations (3.8), (3.11) and (3.12), and let T{ree{x — y), 
T^"^{x — y) andC^\^{x — y) be the corresponding infinite-volume quantities, defined in 
equations (3.18), (3.19) and (3.20). Then the infinite-volume limits of Ti^ee,!v{x,y), 
^w'ir,tv{x^y) and Cy,ir,A{x,y) exist, and 

T^rr(^-y)= lim rfi- (x,y) (3.24) 



and 



C„ir{x -y)= lim Cwir,A(a:, y) . (3.25) 



Remark. For local observables, the existence of the thermodynamic limit follows 
immediately from the FKG monotonicity properties (2.18) and (2.19) — see equa- 
tions (2.20) and (2.21). This, however, does not imply the relations (3.23), (3.24) 
and (3.25), since the events in question are nonlocal; the relations can only be es- 
tablished after an interchange of limits. In the ordered phase, this interchange is not 
merely technical — it is related to the question of how the infinite cluster emerges 
from large clusters in a finite volume. Thus it depends sensitively on boundary con- 
ditions. For example, for a free boundary condition analogue of the finite-cluster 
connectivity (3.19), an infinite-volume statement like (3.24) is actually false. 

Proof. Introducing the event Rx,y{^) that x and y are connected in -B(A), the right 
hand side of (3.23) is 

lim^ /Ufree,A(2^ ^ V) = llHl^ /Ufree,A(^x,y (A)) 

while the left hand side is 

yUfree(a: ^ y) = Hm ;Ufree(^i,y(A')) = ^^^ 1™ /Ufree,A(^x,y (A')) . 

A'^Zd A'^Z^A^Zd 
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We therefore have to show that 

lim Urn ;Ufree,A(^x,y(A')) = hm ;Ufree,A(^x,y(A))- (3.26) 

In order to prove (3.26), we combine the monotonicity property (2.18) with the fact 
that i?x,y(A') C Rx,y{^) if A' C A to get 

/Ufree,A'(i?x,y(A')) < /Ufree,A(i?x,y (A')) < /Ufree,A(i?x,y (A)) if A' C A . (3.27) 

Taking the hniits A — > Z'^ and A' ^ Z'^, the inequahty (3.27) imphes equation (3.26). 
In order to prove (3.24), we consider the event i?^^(A) that x and y are connected 
by a cluster C{x) C -B(A), as introduced in Proposition 2.7. Recalling that dA = 
{x ^ A|dist(a:, A) = 1}, we see that i?^^(A) is the intersection of the event Rx,y{A.) 
with the event that x is not connected to dA. We claim that 

//wir,A'(i?I;"^(A'))<;Uwir,A(i?I;"^(A'))<;Uwir,A(i?^:;,(A)) if A' C A . (3.28) 

As before, the second inequality follows from the fact that R^^y{A') C i?^^(A) if 
A' C A, which implies that /Uwir,A(-Rx^(A')) is monotone increasing in A'. However, 

the monotonicity of /Uwir,A(-Rx"y(A')) ^^ A is less obvious because i?^^(A') is neither 
an increasing nor a decreasing event. It is, however, an event of the form (2.27) 
considered in Proposition 2.6 and its corollary. Namely, 

R%{A') = \J{C{x) = B} = [jiooB = 1} n {uo^-B = 0} , (3.29) 

B B 

where the union goes over all connected sets B C -B(A') that join x to y, ujb is the 
configuration u restricted to the set B, and d*B is the set of all bonds in B^ \ S 
which are connected to B, as in the proof of Proposition 2.7. Thus by the corollary 
to Proposition 2.6, /Uwir,A(-Rx"y(A')) is an increasing function of A' C A, which is 
actually stronger than the first inequality of (3.28). This completes the proof of 
(3.24). 

In order to prove (3.25), we remark that it has been already shown in [ACCN] 
(Theorem 2.3c) that fi^iT-{\C{x)\ — oo) = lim^^zd A*wir,A(3^ ^^ dA) . The proof of 
(3.25) therefore reduces to showing 

A*wir(|C'(3^)| = OO and |C(y)| = oo)) = lim Hwii Aix ^^ dA and y <-* dA) . (3.30) 

A^Zd 

Proceeding as before, we now introduce R^^ as the event that both x and y are 
connected to dA. With this notation, equation (3.30) can be rewritten as 

lim lim ;Uwir,A(^ft ) = lim ;Uwir,A(^f'l) ■ (3.31) 

A'^Z^A^Zd '^ A^Zd '^ 

Using (2.19) instead of (2.18), and observing that R^^ D R^^ if A' C A, we obtain 

^wir,A'(-Rx,y ) ^ ^wir,A(-Rx,y ) ^ ^wir,A(-Rx,y) ^^ A' C A . (3.32) 

As before, the proof is completed by taking the limits A — * Z*^ and A' -^ Z^. D 
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3.2. The Covariance Matrix and its Eigenvalues. 

Here we analyze the structure of the covariance matrix 

Gr{x -y) = {qd{a,, m) ; qd{ay, n)), , (3.33) 

with free and constant boundary conditions, summarizing our results in Theorem 
3.5 at the end of the section. Before discussing particular boundary conditions, we 
note that in general 

J2 Gn^ -y) = j2Gri^-y) = o, (3.34) 

m n 

which follows from the fact that any truncated expectation {A; S)^ vanishes if either 
A or S is constant, and from the obvious relation X^meS ^i'^x^'m) = 1. In particular, 
this implies that, independent of boundary conditions, G]^^ always has a trivial 
eigenvalue 0, corresponding to an eigenvector vq = {!,■ ■ ■ ,1) E M^. 
Now consider the matrix with free boundary conditions 

G^ei^-y) = {qS{a.,m);qd{ay,n)),^^^. (3.35) 

Due to the permutation symmetry of the Hamilton function (2.1) and the symmetry 
of the boundary conditions, all diagonal elements are equal, as are all off-diagonal 
elements. Combining this with the observation (3.34), we conclude that 

G^elix -y)= (^S{m, n) - (1 - 6{m, ^))^) G'fK^ - y) ■ (3-36) 

Given (3.36), the matrix G'^^^{x — y) is easily diagonalized. We find one trivial 
eigenvalue 0, corresponding to an eigenvector vq = (1, ■ ■ ■ ,1), and one {q — l)-fold 
degenerate eigenvalue 

Gi,ee{x -y) = — yC^reela^ ' V) = — ^(^^(^X' ^v) " l)free ' (3-37) 

corresponding to the {q — 1) -dimensional eigenspace orthogonal to vq. In the second 
equality in (3.37), we have reexpressed G^^^^{x — y) as the usual two-point function. 

Remark. The above results imply that, in the free boundary condition case, the 
covariance matrix of the g-state Potts model does not contain more information 
than the standard two-point function. As we will see below (and as should be 
clear from the fact that G'^'^{x — y) always has one trivial eigenvalue), the same is 
true of the covariance matrix of the Ising model (q = 2) with constant boundary 
conditions. This may explain why the covariance matrix has not been more widely 
studied previously. However, as we shall see below and in subsequent sections, the 
q > 3 state matrix with constant boundary conditions does have additional content, 
and this content has a clear stochastic geometric interpretation. 
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Next we analyze the covariance matrix with constant boundary conditions, 

Gr{x -y) = {q6{a^, m) ; q6{ay, n))^ . (3.38) 

Starting with the special case g = 2, we use (3.34) to conclude that G^J^{x — y) = 
Gl^{x-y) = -C^^^ix-y) = -Gl^{x~y). Combined with the fact that Gl]" (a; -y) = 
G\^{x — y) by the symmetry of the model, we obtain 

G]^''{x-y) = {d{m,n)-{l-6{n,m))G'^o^{x-y) forq = 2. (3.39) 

Observing that the matrix structure of (3.39) is identical to that of (3.36) with 
q = 2, we see that we again obtain a trivial eigenvalue of and an eigenvalue 

<^wil(^-?/) = G^o°(^-2/) = («x;Sy)o forQ = 2. (3.40) 

Here we have rewritten GQ*^(x — y) in terms of standard Ising spins s^ = 2d{ax,0) — l. 
For q ^ 2, the matrix structure of G^"^(x — y) is less trivial. Using relation 
(3.34) and the fact that constant boundary conditions c & S leave the symmetry 
of permutations among elements of 5" \ {c} unbroken, it is easy to show that there 
are only two independent matrix elements. Taking these to be Gq^{x — y) and 
Gq^{x — y), we obtain 

G''f{x-y) = G'^o^{x-y) if ceS, 
Gr{x-y) = Gl\x-y) if n^c, 

G':^(x-y) = G';^^{x-y) = —G9.^(x-y) if n^c, and 

q-1 

Grix - ^) - (g _ i|(g _ 2) ^o°(^ -y)- ^Gl\x - y) if n^m^n^m^c. 

(3.41) 

In order to diagonalize G'^^{x — y), we begin by observing that the expectation 
(■)^ is invariant under the group Sq-i of permutations of 5'\{c}. Diagonalizing G'^^ 
first on the Hilbert space corresponding to the trivial representation of Sq-i, we 
identify two eigenvectors: vq = {1, ■ ■ ■ ,1), corresponding to the simple eigenvalue 
zero, and Vi, with components (i'i)m = q5{m, c) — 1, corresponding to the nontrivial 
simple eigenvalue 

G^2{x -y) = -;^G'o\x ~y) = -^{q5{a,, 0); q5{ay, 0)), . (3.42) 

On the remaining {q — 2) -dimensional subspace orthogonal to vq and vi, we finally 
obtain the {q — 2)-fold degenerate eigenvalue 

G^ix -y) = ^^ Gl\x -y)- ^^ _ ^^^^ _ ^^ Gl\x - y) 

= Gl\x-y)-G]^{x-y). (3.43) 
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It is interesting to note that, as in (3.37), it is possible to express the eigenvalue 

(2) 

G^jj.(a: — y) in terms of an nntruncated expectation. Indeed, we may simply rewrite 
the second line in (3.43) as 

gS(x -y) = liiqSia,, 1) - q6{a,, 2)) ; {q6{ay, 1) - q6{ay, 2))), 

= \{iq6ia^, 1) - q6ia^, 2)) {q6{ay, 1) - q5((7„ 2)))^ . 

^ (3.44) 

We both summarize the results of this section and establish their stochastic 
geometric significance in the following: 

Theorem 3.5. Consider the q-state Potts model with q > 2. Then the free bound- 
ary condition covariance matrix G^^^{x — y) has the simple eigenvalue zero corre- 
sponding to the eigenvector vq = (1, ■ ■ ■ , 1), and a (q— l)-fold degenerate eigenvalue 

Gireeix - y) = — y (?5((T^, (T^^) - 1)^^^^ (3.45) 

corresponding to the subspace orthogonal to vq . For q > 3 and all c E S, the constant 
boundary condition covariance matrix G'^^{x — y) has the simple eigenvalue zero 
corresponding to the eigenvector vq = (1, ■ ■ ■ ,1), a nontrivial simple eigenvalue 

G^2i^-y) = -^{qS{a,,0);qS{ayMo (3-46) 

corresponding to the eigenvector vi, with components {vi)m = qS{m,c) — 1, where 
Vq and vi belong to the trivial representation of the unbroken subgroup Sq-i, and a 
(q — 2)-fold degenerate eigenvalue 

Gwil(^ -y) = ^((Q'^K, 1) - ?5(a,, 2)) {q5{ay, 1) - q5{ay, 2)))^ (3.47) 

corresponding to the subspace orthogonal to vq and vi. For q = 2, the matrix 
Qrnn^^ — y) has Only the trivial eigenvalue zero and the eigenvalue G\^-Ij.{x — y). 

Moreover, the eigenvalues Gf^^eix — y), G\^(^{x — y) and G^(^{x — y) can be ex- 
pressed in the random cluster representation as 

Gfreeix -y) = qnree{x-y) , (3.48) 

(^will^ -y) = Qr^Ux -y) + q{q - l)C^ir{x - y) (3.49) 

and 

G^2{^-y) = qr^Xx-y). (3.50) 

Proof. It only remains to establish the random cluster representations (3.48), (3.49) 
and (3.50) of the eigenvalues. But these follow immediately from expressions (3.37), 
(3.42) and (3.43) for the eigenvalues in terms of the matrix elements, and expressions 
(3.21) and (3.22) relating the matrix elements to the random cluster connectivities 
and cluster covariance. D 
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4. The Correlation Lengths 

4.1. Existence of the Lengths ^^ee, C^l and ^^f^. 

In this subsection, we establish the existence of the limits (1.7), (1.8) and (1.9) 
using standard reflection positivity arguments. Namely, introducing the unit lattice 
vector ei = (1, 0, ■ ■ ■ , 0) G Z'^, we prove the following: 

Theorem 4.1. Let q > 2 be an integer, and let G{t) denote G{ree{tei), G^/j.(tei) 

(2) 

or (for q > 3) G)^-^^{tei), see Theorem 3.5. Then G{t) is a non-negative, mono- 
tone decreasing and log convex function oft, so that {G{t)/G{0))^^^ is monotone 
increasing in t, the limit 

l/?--limi^ (4.1) 

exists and the function G{t) obeys the a priori bound 

G{t) < G(0)e-*/«. (4.2) 

Furthermore, denoting the correlation lengths defined in (4-1) by ^iree, ^wir ^'^^^wir' 
we have 

e^^.^ > e^^\ (4.3) 

Proof. We start with the observation that each G{t) can be written as a truncated 
(inflnite- volume) expectation of the form 

G{t) = {A;T'A),, (4.4) 

where b denotes either free or constant boundary conditions, A = A{ax) is an 
observable which depends on the spin variable a^ of a single point x E Z^, and T^A 
is the translation of the observable A by tei. Equation (4.4) follows from (3.37) and 

(3.35) for Gfree, from (3.42) for G^V and from (3.44) for G^l. 

Due to the reflection positivity of the model (see Appendix A for a review of the 
basic ideas), T can be represented as a non-negative contraction (0 < T < 1) on a 
Hilbert space H, and 

G'(t) = (^,TV) 

for a suitable vector tf) E Ti. Obviously, this implies that G{t) is a monotone 
decreasing, non-negative function of t. By the Cauchy-Schwarz inequality, 

G{\{t,+t2)) = (T*i/V,T*^/'V') < VG{ti)G{t2), 

so that G{t) is log convex. Noting < G{0) < oo, this implies that (G'(t)/G'(0))^/* 
is monotone increasing in t, which in turn immediately implies existence of the limit 
and the a priori bound. Finally, (4.3) follows immediately from the representation 
in Theorem 3.5 and existence of the limits. D 
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Remark. For Gfree and G^/j.' ^^^ existence of the corresponding correlation lengths 
can also be established by subadditivity arguments (see Section 4.2 below). While 
these arguments are more involved than the reflection positivity proof presented 
above, they have the advantage that they give existence of limits analogous to (1.7) 
and (1.9) for non-integer values of g, defined directly in terms of Tfree and r^"^. 
(c.f. (equations (3.48) and (3.50)). Moreover, a slight variation of these arguments 
can be used to establish left continuity of the inverse correlation length l/^free(/3) 

(2) 

(Theorem 4.2) and upper semicontinuity of l/^^^ir (Theorem 4.3). On the other 
hand, subadditivity does not establish log convexity of G{t) and hence monotonicity 
of the full sequence (G(t)/G'(0))^/*, as we have from the above theorem for integer 
q. Furthermore, we are not aware of any proof of the existence of ^^j^. which does 
not involve reflection positivity. 

4.2. Equivalent Characterizations of ^free? ^wir ^^^ Wir* 

We already have stochastic geometric representations for the correlation lengths 
^free and ^^jj. as the decay rates of of r^ee and r^"^. - see Theorem 3.5. In this 
subsection, we provide a stochastic geometric representation for ^^j^. (Theorem 4.4) 

(2) 

and give alternative representations for ^free and ^^j^. (Theorem 4.3, Lemmas 4.6, 
4.7 and 4.8). On the one hand, these alternative representations allow us to prove 

(2) 

several results on the behaviour of ^free and ^^^^, in particular left continuity of 

(2) 

l/^free(/3) (Theorem 4.2), upper semicontinuity of l/^^jj. (Theorem 4.3), and the 

(2) 

two-dimensional dichotomy (1.11) and (1.12) involving ^free and ^^j^. discussed in 
the introduction. On the other hand, the alternative representations may be of 

(2) 

interest for numerical determinations - in particular the representation of ^^j^. in 
terms of the probability T^\'^^{n) that the diameter of the cluster C(0) is n (The- 
orem 4.3). The representation of ^^jj. as the decay rate of the covariance C^ir, 
provided the magnetization M(/3) > 0, may be of interest both to mathematicians 
and numerical physicists. It is worth noting that many of the results of this subsec- 
tion are generalizations of corresponding percolation results of [CCGKS] to q > 1, 
but the proofs are quite different due to the lack of independence, the lack of a BK 
inequality, and the presence of boundary conditions. 

Some of the results in this section (and most of the proofs) are of a rather techni- 
cal nature. In particular, we introduce many connectivity functions and ultimately 
show that they have only a few independent decay rates. However, in the process, 
the notation and the arguments become rather cumbersome. In order to simplify 
matters, we flrst introduce only a few "physical" connectivity functions and sum- 
marize the results of independent interest on ^free, ^wir ^^^ Wir i^ Theorems 4.2, 4.3 
and 4.4, respectively. The remainder of the subsection is devoted both to the proof 
of these results and to the statement and proof of several more techncial results 
which we will need for our proof of the two-dimensional dichotomy in Section 5. 

We start with a few deflnitions. For b = "wir" or "free" , we introduce the on-axis 
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connectivity function 

Tb{Lei)=iJb{0^ Lei), (4.5) 

the on-axis finite-cluster connectivity 

ri^iLei) = fib{{0 ^ Lei} n {|C(0)| < oo}), (4.6) 

the diameter function 

T^'^'^iL) = ;U,(diamC(0) = L), (4.7) 

where diamC(O) denotes the diameter of the cluster C(0) in the i^o norm, i. e. the 
maximum diameter in any of the d coordinate directions, and the covariance 

a{x - y) = fiB{{\C{x)\ = oo} n {\C{y)\ = oo}) - PliPf, (4.8) 

where 

P^(/?) = ^5(|C(0)| = oo). (4.9) 

We denote the corresponding correlation lengths - whenever they exist - by ^b, ^f ", 

^diam g^^^ ^C_ 

Theorem 4.2. Let < /3 < oo and q > I. Then the correlation lenths ^wir o,nd 
^free exist, ^free < ^wir; o,nd l/^free is a left continuous function of 13. 

Theorem 4.3. Let < (3 < oo and q > I. Then for b = "wir" or "free", the 
correlation lengths ^^" and ^^^^^ exist and are equal: ^f^ = ^^i^™. Also ^^? < ^fj-gg, 
and 1/^wh- ^-^ ^'^ upper semicontinuous function of /3. If q > 3 is an integer, then 
in addition 

t(2) ^ cfin ^ ediam (a -i r)\ 

^wir Swrir S«rir • V "/ 



Remark. Combined with the obvious inequality ^f^^ < ^free, the inequalities from 
Theorem 4.2 and Theorem 4.3 give 



4 wir — Sfree — sfree ^ 4 wir 

provided < (3 < oo and q > I. 

Theorem 4.4. Let < (3 < oo and q > 2 be an integer. Then the correlation 
length ^^^^ exists and 

eir = Cl ^f M(/3)>0 (4.11) 

while 



Cr = and ei^.J = e^2.J = ^free tf M{(3)=0. (4.12) 



In order to prove Theorem 4.4, we use a proposition which may be of independent 
interest and is stated next. 
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Proposition 4.5. Let < (3 < oo and q> 1. Then 

a{x-y)>ri-{x-y)PUP). (4.13) 

Proof of Proposition 4-5. The proof of this proposition is an easy generalization of 
the proof of the corresponding statement in [CCGKS]. For a set S C Md, let P{B) 
be the set of points x such that x E db for some bond b E B. Denoting by B^ the 
family of finite connected subsets B cM^ for which x G P{B), we have 

Ct{x -y)= ^i,{\C{x)\ < oo)^,(|C(y)| = oo) - ^Ji,{{\C{x)\ < oo)} n {\C{y)\ = oo}) 
= Y. if^biCix) = B)^,{\C{y)\ = oo) - M{C{x) = B} n {\C{y)\ = oo})) 

= Y^ ^biC{x) = B){^,{\C{y)\ = oo) - M\C{y)\ = oo | C{x) = B)) 
BeB^ 

> J2 MC{x)=B){^t{\C{y)\ = oo)-M\C{y)\ = oo\C{x)=B)) 

y(=P(B) 

= J2 f^b{C{x)=B)fi,{\C{y)\ = oo)=rt{x-y)P'^{(3) 



BeBa;: 

yeP(B) 

where in the fourth step we have used that for all B E Bx 

{C{x) =B} = {UB = 1} n {ud^B =0} = A2nD 

is an event of the form considered in the second inequality (part 2i) of Proposition 
2.6, and hence 

^ib{\C{y)\ = oo) - ^ibi\Ciy)\ = oo | C(x) = B)>0. 

Here, as in the proof of Proposition 2.7, ujb is the configuration uj restricted to B, 
and the boundary d*B of B is the set of all bonds in B^ \ i? which are connected 
to B. D 

Proof of Theorem 4.4. Since M(/3) = PZ''W), we have that for M(/3) > 
{q - l)Cwir(x -y)+ T^Ux -y)< ({q - 1) + ]^^)Cwir(a; - y) < 

< (1 + (^— 1]mm) ^^^ " ^^^""^"^ " ^^ ^ ^'^^^"^ " ^^^ 

by Proposition 4.5 and the fact that r^^ > 0. Combined with Theorem 4.1., which 
guarantees the existence of the inverse correlation length 

. /AD ^ _ ,. logjiq - l)C^ULei) + r^ULei)) 

/ ?wir llUi 

L — >oo J_/ 
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provided g > 2 is an integer, we obtain the statement of Theorem 4.4 for M(/3) > 0. 
On the other hand, if M{/3) = P^'{I3) = 0, then t^I = Twir and C^ir{x - y) = 
which imphes ^^^l = ^^^l and ^^^^ = 0. FinaUy, M{(3) = imphes //wir = yWfree 
[ACCN] and hence ^gj = ^ree- □ 

In order to prove Theorems 4.2 and 4.3, we will need several approximations to 
the connectivity functions Tb{Lei) and r^'^(Lei). Additional approximations will 
be needed to prove the dichotomy (1.11) and (1.12) discussed in the introduction. 
Rather than introducing them as they arise, we define all of them here, so that 
the reader may more easily refer back to the definitions. We will consider several 
subsets of Z'^, namely the the "cylinder" 

H{L) = {xeZ'^\0<xi<L}, (4.14) 

the "tunnel" 

T{M) = {xeZ'^\-M/2<x^<{M+l)/2 i = 2,---,d} (4.15) 

and the "box" 

A{L,M)=T{M)nH{L). (4.16) 

We then consider the following approximation to TbiLei) in the cylinder, tunnel 
and box: 

r,^^^(L)=^,,H(L)(0--Lei), (4.17) 

rl'^^'iL^M) = fib,nM){0 ^ Lei) (4.18) 

and 

r^^'^'iL^M) = ;Ub,A(L,M)(0 ^ Lei). (4.19) 

In all cases b = "wir" or "free". Assuming that they exist, we denote the corre- 
sponding correlation lengths by ^^^ , ^^""(M), and ^^°^(M). We also consider the 
several approximations to r^"(Lei), namely 

f',^\L) = ;U,,H(L)({0 ^ Lei} n {|C(0)| < oo} n {0 /. dH{L)}), (4.20) 

f^-(L, M) = ^t,A{L,M)i{0 ^Lei}n{0^ c>A(L, M)}), (4.21) 

f;;^\L) = iib{{0 ^ Lei} n {|C(0)| < cx)} n {C(0) C B{H{L))}) (4.22) 



and 



'b 

ibox 



fr"(L, M) = ^,({0 ^ Lei} n {C(0) C B{A{L, M))}), (4.23) 

denoting the corresponding correlation lengths - whenever they exist - by ^^^ , 
4box(M),ff and4t-ox(^). 

We note that the distinction between (4.17) - (4.21) and (4.22), (4.23) is that in 
the former quantities the probabilities are computed with respect to measures that 
live on the relevant sets A, while in the latter the probabilites are computed with 
respect to the full measures fib, but the events in question occur in the relevant sets 
B{A). 

Our first lemma gives the equivalence of several definitions of the correlation 
length ^b and will be used at the end of this section to prove Theorem 4.2. 
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Lemma 4.6. Let < (3 < oo and q > I. Let t{L) denote T^ir{Lei), rfree(-^ei), 
rj*."g(L, M), Tf*j.gg(L), or r^°^{L, M). Then the correlation length ^ corresponding to 
t{L) exists, and t{L) < e~^'^ . Furthermore, the correlation lengths ^\^^e{M) and 
^}?°g(M) are monotone decreasing in M , ^f^el-^) = ^free(-^)' ^'^^ 

^free = ^f^e = ^free = ^free^ (4.24) 

where ^7, = limM^oo Ce(^) <^nd HZ = liniM^oo ilTe{M). 

Proof. Considering an arbitrary subset A C Z^ and two points x and y in A, we 
note that by the FKG inequahty (2.17) 

lJ'h,K{x ^y)> fJ^bM^ ^ z)iJb,Ai^ ^ y) (4.25) 

for all 2 G A; furthermore, by the FKG monotonicity (2.18) 

yUfree,A(a^ "-^ V) > IJ'hee,A'{x ^ y) (4.26) 

for all A' C A containing x and y. Using these inequalities, one obtains subaddi- 
tivity, and hence existence of the corresponding correlation length ^, together with 
the a priori bound t{L) < e~^'^ for all five connectivity functions t{L) considered 
in the theorem. Observing that the monotonicity (4.26) implies the monotonicity 
of Tf^ee (^' ^) and tI^^^{L, M) in M, one obtains the monotonicity of i\^^^{M) and 
^}?°g(M) in M, as well as the justification of the interchange of limits 

M^oo L^oo L L^oo M^oo L 

and similarly for t^°^- The only additional ingredient needed in the proof of the 
equalities 

^free = lim Ce(M) and Clt = Jini ClTe(M). 

is that A^free, T(M) (3^ ^^ v) couverges to liireei^ *-^ y) (and similarly for t^°^), which 
is established in the same way as (3.23). 

We are left with the proof of the equalities ^free = CfrL and ^f"ee(^) = ^free(^)- 
To this end, we use (4.25) and (4.26) to get the bound 

n-l 
/Ufree(0 ^ uLei) > /Ufree,H(nL) (0 ^ TlLci) > J]^ /Ufree,H(nL) («-^ei ^ (z + l)Lei). 

i=0 

(4.27) 
Taking the limit n ^ oo, and noting that all but say ^/n terms on the right hand 
side have arguments which are sufficiently far from the boundary, we obtain 

g-i/Cfroo > ^-L/^t^l > rfree(^ei) 
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which, in the hmit L -^ oo, imphes that ^free = ^free- '^^^ equahty of ^fj!^e(Af ) and 
^f^°g(M) is proved in the same way. D 

The next two lemmas give several useful relations between the correlation lengths 
corresponding to (4.20) - (4.23), and are important ingredients for the proof of 
Theorem 4.3 (see below) and for the proof of the dichotomy (1.11) and (1.12) (see 
Section 5). In order to state the first of these two lemmas, we introduce for each 
x±_ G Z^~^ n [— M, MY~'^ the off-axis connectivity function in the box, 

f^^(L, M; x^) = ^5({0 ^ (L, x^)} n {C(0) C S(A(L, M))}), (4.28) 

and for each x_l G Z'^"^ the off-axis connectivity function in the cylinder 

ff{L-x^) = in,{{Q ^ (L,Xi)} n {|C(0)| < oo} n {C(0) C B{H{L))}). (4.29) 

We note that Am = {C'(O) C B[K{L, M))} is an increasing sequence of events which 
converges to the event {|C(0)| < 00} fl {C(0) C B{H{L))}. As a consequence, 

fb^°^(L, M; x±) / f^^\L; x±) as M / cx). (4.30) 

Lemma 4.7. Let < (3 < 00 and q > I. Then for b = "wir" or "free", the 
correlation lengths $,^^ and ^j^°^(M) corresponding to the connectivity functions 
(4.22) and (4.23), as well as the limit l^""" = limM^oo lb°''(M), exist and 

if = 1^^ (4.31) 

In addition, 

fr(L,M;x^) < C(/3,Q)e-^/«"°''W < C(/3, Q)e-^/^"°^ (4.32) 

and 

ff{L-x^)<C{fi,q)e-'^'^l'' (4.33) 

where C{(3,q) < 00 is continous as a function of 13 and independent of L and M. 

Remark. We will later show that ^j;"^^ = ^^" (see equation (4.53)), so that Lemma 
4.7, together with Theorem 4.3, gives us yet another characterization of ^4^;^.- 

Proof. In order to prove the lemma, we first establish the existence of a constant 
C(/3, q) <oo such that for all x^ G Z'^'^ n [-M, M]^"!, 

f^^(Li, M; x,_) f^"(L2, M; x,_) < C(/3, q) f^'^(Li + L2 + 1, M; 0). (4.34) 

To this end, we first rewrite the left hand side of (4.34) in a form which allows us 
to apply Proposition 2.7. We consider the boxes Ai = {x G T{M) | < xi < Li}, 
A2 = {a; G T{M) | Li + 1 < xi < Li + L2 + 1} and A = A(Li + L2 + 1, M), and 
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the points X- = {Li,x±)j Xj^ = (Li + l,x±) and y = (Li + L2 + 1,0). Using the 
translation and reflection invariance of the measures Hb, we then rewrite 

f^^(Li, M; x^) f^^(L2, M; x^) = ^ib{R^%_ (Ai)) fib{R^l^yiA2)) (4.35) 

where, as in Proposition 2.7, i?o"_(Ai) is the event that and X- are connected 
by a cluster C(0) C S(Ai), and similary for R^'^y{A2). Observing that S+(Ai) n 
B{A2) = B{Ai) n S+(A2) = 0, we apply Proposition 2.7 to obtain 

/x,(i?fi';_(Ai))^,(i?fi-,^(A2))<^,(i?|}",_(Ai)ni?fi^,^(A2)). (4.36) 

Next we note that all configurations u E i?o'^_(Ai) fl -Rx+,y(^2) would contribute 
to the event RQ^y{A) if the vacant bond (a:_, x_|_) were occupied. Using finite energy 
in the form (2.24), we therefore conclude that 

^^biRt_{A^)nR'^UA2)) < C{(3,q)^ib{Rty{A)) (4.37) 

for a suitable constant C(/9, q) < 00. Observing that 

^ib{Rty{A)) ^ f^°^{L, +L2 + I, M; 0) , 

we obtain the subadditivity bound (4.34). By standard arguments, the bound 
(4.34), together with the monotone convergence (4.30), implies the existence of the 
correlation lengths i'^^ and fj;'°^(M), and the hmit fj;'"^ = hm^^oo ^^"^(Af ), the a 
priori bounds (4.32) and (4.33), and the equality of ^l^ and ^b°^. Finally, we note 
that by the finite energy relation (2.24), without loss of generality C(/9, q) may be 
chosen to be a continuous function of (3. D 

Lemma 4.8. Let < /3 < 00 and q > 1. Then the correlation lengths ^^^^ and 
$,^°^{M) corresponding to the connectivity functions (4-20) and (4-21), as well as 
the limit ^^"f = limM^oo ^w?:^(^)j exist and 



In addition, 



and 



S = ^ = i^- (4.38) 



r^r(^, M) < C{/3, Q)e-^/^w°?(M) < c(p^ q)e-^/^^ (4.39) 



r:^;(iv)<C(/3,Q)e-^/C^^ (4.40) 

where C{(3j q) < 00 is continuous as a function of (3 and independent of L and M . 

Remark. While Lemma 4.7 gives us alternative representations of ^^" in terms of 
decay rates of infinite-volume connectivities, this lemma - together with equation 
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(2) 

(4.53) and Theorem 4.3 - gives us representations of ^^j^. in terms of finite-volume 
quantities. 

Proof. Following the proof of Lemma 4.6, we first establish two inequalities analo- 
gous to (4.25) and (4.26). In order to state them, we introduce the cylinders 

if (Li, L2) = {x G Z^ I Li < xi < L2}, (4.41) 

the boxes 

A(Li, L2, M) = H{L^, L^) n T{M), (4.42) 

(with T{M) as defined in (4.15)), the events 

Rll^A) = {^1^1 ^ ^2ei} U {C{0) C S(A)} , (4.43) 

and in particular 

RIIl,{M) = RllLMiLi,L2,M)). (4.44) 

We then claim that for a suitable constant C{(3, q) < 00, 

^^^ir,AiRllL,iM))^^^.r,AiRt+l,L,iM)) < C {^ , q) ^^,,,K{Rt,L,{M)) (4.45) 

if AdA(Li,L2,M), 

/^wir,A'(i?S,L.(A)) < ^wir,A"(i?i'J,L.(A)) (4-46) 

if A" D A' D A, and 

^wir,A(i?S,L.(A')) < ^wir,A(i?i^,L.(A")) (4-47) 

if A' C A" C A. While the first of these three inequalities is proved in exactly the 
same way as (4.34) using finite energy and Proposition 2.6, the second follows from 
Proposition 2.6 (and its corollary) alone since i?£" ^ (A) is an event of the form 
(2.27), see proof of Proposition 2.7 and equation (3.29). The last of the inequalities 
follows from the fact that Rf^j^^iA') C R11l^{A") if A' C A". 

Given (4.45) - (4.47), the proof of Lemma 4.8 is analogous to that of Lemma 4.6, 
with 

C(/3,Q)^-Vwir(^0,(nL+n-l)(M)) 

> C{l3,qY fJ'wiT,A{0,nL+n-l,M){Ro,{nL+n-l)iM)) 
n-1 

> 11 fJ'wir,A{0,nL+n-l,M)(Ri{L+l),L+i{L+l){M)) (4.48) 
i=0 

replacing the inequality (4.27). D 
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We finally turn to the proofs of Theorems 4.2 and 4.3. 

Proof of Theorem 4-2. The existence of the correlation lengths ^free and ^wir has 
already been established in Lemma 4.6, and the inequality ^free < ^wir follows 
immediately from the FKG ordering (2.23), so all that remains to show is left 
continuity of I / CireeiP) ■ Due to equation (4.24), l/^free(/3) is a limit of finite- volume 
(and hence continuous) quantities, namely 

!/«,„,(« ^ - lim ,to l^i:tf^i£). (4.49) 

M — >oo L — >oo Li 

As shown in the proof of Lemma 4.6, the finite-volume connectivity t^°^{L, M) is 
subadditive in L and monotone increasing in M. It is also monotone increasing and 
continuous in (3. Choosing suitable subsequences, e. g. L = 2", and noting the minus 
sign in (4.49), this gives l/^free(/3) as the limit of a decreasing sequence of continuous 
decreasing functions, and hence establishes the desired left continuity. D 

Proof of Theorem 4-3. We start with the obvious bounds 

f,^-(L, M) < r/\L) < rt{L) < J2 rt'^in) (4-50) 

n>L 

and 

f^^°'^(L, M) < rf ^"(L) for all M < L. (4.51) 

Now consider the event {C(0) = B} where B is some given set of diameter L. Using 
suitable rotations and translations by vectors t G Z'^, |t| < L, each such cluster B 
can be transformed into a cluster B C A(L, 2L) connecting the origin to a point 
X = (L, x±) in the boundary of A(L, 2L). As a consequence. 



.diam(^) < d{2L + 1)^ V f6^°"(L, 2L, x^) < d{2L + i)^d-i^-L/i,- ^ (452) 






where we have used the a priori bound (4.32) of Lemma 4.7 in the last step. 

Combining the bounds (4.51) and (4.52) with Lemma 4.7, we immediately obtain 
the existence of the correlation length ^^lam ^^^^ ^y^^ equality of ^^'^™ and ^5°^. 
Combining the bounds (4.50) and (4.52) gives the existence of ^^" and the equality 
of ^^" and ^^"^, provided ^^"^ < 00. If, on the other hand, ^}l"^ = ^|j^ = 00, we use 
the bound f^^ (L) < t^'^{L) < 1 to prove that the inverse correlation length 1/^f" 
exists and is equal to zero. Thus we have the existence of the correlation lengths 
^?" and ^u'^'^^ and the equality 



3m 4.3, namely ^^^ 
immediately from relation (3.50) of Theorem 3.5. 



f2) ft 

The final equivalence of Theorem 4.3, namely ^^^^ = ^^^^ for integer g > 3, follows 
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We are therefore left with the proofs of the inequahty ^^P < ^f^^ and the upper 
semicontinuity of l/^^^j.- Noting that r^" is the probabihty of an event of the form 
considered in Proposition 2.6, the inequahty follows immediately from the infinite- 
volume limit of (2.31). To prove the upper semicontinuity, we note that by Lemma 
4.8 and equation (4.53) above, l/^^h- '^^^ be written as a limit of finite-volume 
quantities, namely 

VC^r = yi^r = yC^^ = - Jim lim ^ --/ ' ^ . (4.54) 

Combined with the a priori bound (4.39) of Lemma 4.8, this implies 



e '^""^ = sup 






and hence 



l/e=M{'°«^<^'«'-^°«"-'^'^^'}. (4.65) 



Since both logC(/3, g) and — log'r^?'^(L, M) are continuous and hence upper semi- 
continuous functions of /3, and since the infimum of upper semicontinuous functions 
is upper semicontinuous, this establishes the upper semicontinuity of l/^^^r- ^ 

In the above proof, we actually obtained one additional equivalence which is not 
stated in Theorem 4.3, but which will be necessary in the proof of our dichotomy. 
Namely, by equation (4.53), we have: 

Corollary. Let < P < oo and q > 1. Then 1^°^^ = ^^". 

5. The Two-Dimensional Dichotomy and Related Results 

5.1. Heuristics and Preliminaries. 

The goal of this section is a proof of the two-dimensional dichotomy, the principal 
part of which is the duality relation (1.11) for all /3 in the low-temperature regime. 
In this subsection, we discuss the heuristics of the relation, state our results and 
briefiy review two-dimensional duality in the random cluster model. In the next two 
subsections, we derive upper and lower bounds on r^P and its approximations in 
terms of Tfree and approximations to r^^^- Finally in the fourth subsection, we put 
these bounds together with the equivalence lemmas of Section 4 and the ergodicity 
theorem of Section 2 to prove the dichotomy. 

In order to explain the heuristics of the duality relation (1.11), let us consider the 
representation of the random cluster model in terms of the order-disorder contours 
introduced in [LMMRS], see also [BKM]. In this representation, contours are defined 
as (the connected components of) the boundaries between regions of occupied bonds, 
regarded as ordered regions, and those of vacant bonds, regarded as disordered 
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regions. Notice that in the wired measure, any finite cluster of occupied bonds 
must be separated from the infinite occupied cluster by a (disordered) region of 
vacant bonds. Thus all configurations contributing to T^fj-{x — y) have at least two 
contours surrounding the points x and y - one being the boundary between the 
cluster connecting x and y and the disordered region, and the second being the 
boundary between the disordered region and the infinite cluster. 

Let us begin by considering systems with first-order transitions at the transition 
point [3o- Since both the ordered and disordered phases are stable at [3o, the two 
contours need not remain near each other. Indeed, under similar circumstances, 
[MMRS] proved that two such order-disorder interfaces tend to behave like inde- 
pendent interfaces, leading to a surface tension aoo between two ordered phases 
which is exactly twice the surface tension aod between an ordered and a disordered 
phase. Now in our case, the minimal combined area of the two interfaces is 4|a: — y|. 
Moreover, we expect that large interfaces are suppressed at a first-order transition. 
Thus we expect T^i^(x — y) to decay exponentially with a rate Aaod = '^o'oo, which 
would imply 

l/^^^ = 2aoo. (5.1) 

Obviously, this relation should also be satisfied trivially at (3o for systems with 
second-order transitions - both sides should vanish. 

Now consider the regime /? > /9o in a system with either a first- or second-order 
transition. In this regime, the disordered phase is unstable, so that large regions 
of vacant bonds are supressed. Thus the two contours surrounding x and y tend 
to bind together, leading to a single order-order interface surrounding the points 
of minimal area 2\x — y\- This leads to a exponential decay with a rate 2a"oo, and 
hence again the relation (5.1). 

Note that due to the duality relation aooil3) = l/^free(/3*), equation (5.1) is 
equivalent to the desired relation (1.11). 

It would be interesting to make the above heuristic arguments rigorous. While 
this could presumably be done for sufficiently large q, a direct translation of these 
heursitics into a proof for arbitrary q seems much more difficult. We therefore follow 
a different route, based on our inequalities involving decoupling events (Proposition 
2.6) and the equivalences established in Section 4. 

Before stating our main result, let us recall that the dual inverse temperature (3* 
is defined by 

(e^-l)(e^* -1) = Q. (5.2) 

Our main result is: 

Theorem 5.1. Let d = 2, q > 1 real and < (3 < oo. Then either 

P^-(/3*) = and e5^(/3) = i6ree(/3*) (5.3) 

or 

Pfr^^(r)>0 and CliP) = CireeiP) . (5.4) 
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Remarks. 

1) If, in addition, g > 3 is an integer, it follows easily from the results of the last 
section and the duality relation (1.13) (a proof of which is given in Section 5.4) 
that (5.3) and (5.4) may be replaced by the dichotomy: Either 

P^-(/3*) = and ^iii(3) > e2i(/3) = i6ree(/3*) (5.3') 

or 

Pfr'^^(r)>0 and Cl'lm = £liP) = ^ireeiP) . (5.4') 

2) It follows from the duality relation (1.13) and the monotonicity of P^ {h = "free" 
or "wir") as a function of (3 that the first branch of the dichotomy (i. e. (5.3) or 
(5.3')) occurs when [3 > (3o, the self-dual point, and that the second branch (i. e. 
(5.4) or (5.4')) occurs whenever (3 < (3t = mf{(3 \ M{(3) > 0}. It is presumably 
the case that (3o = /9t, but rigorously this is only known for sufficiently large q 
(see e.g. [LMR] where this has been shown for q > 25). 

We close this subsection with a few remarks on duality. As usual, the dual 
site lattice (Z*)^ is the set of points x* = {xl,X2) G (Z + |)^ with half-integer 
coordinates, and the dual bond lattice 182 is the set of nearest neighbor bonds in 
(Z*)^. To each bond b EM2, there corresponds a dual bond b* G B2 which has the 
same midpoint as b. Similarly, to each configuration w on S C B2, there corresponds 
a dual configuration uj* on B* = {b* \ 6G-B}, given by 

r if uj(b) = 1 , , 

^ (b) = < 5.5 

^ ^ \ 1 if uj{b) = 0. ^ ' 

We will sometimes refer to the bonds b* G B2 for which uj*{b*) = 1 as occupied dual 
bonds. Given a finite box 

A = {x G Z^ I < xi < L, -M/2 <X2<{M + l)/2} 

and the corresponding set of bonds i?"'"(A), one defines the dual of A as 

A* = {x G (Z*)2 I 3y G (Z*)2 with (x, y) G (S+(A))*}. (5.6) 

Note that in general A and A* are not of the same cardinality. Using the appropriate 
Euler relation to relate #(t<;) to #(a;*), it is straightforward to check that for a given 
configuration u on i?+(A) and its dual cu* on i?(A*), 

Gw\t:,I3,k{u}) = G'free,/3*A*('^*) , (5.7) 

where we have explicitly indicated the temperature dependence of the weights (2.14) 
and (2.16). Thus for each A G J-'b+{A)j 

/Uwir,/3,A(^) = /Ufree,/3*A*(^*), (5.8) 

where A* is the event A* = {u* \ u G A}. In the next two subsections, we will often 
characterize events A G J^b+{A) ^^ terms of the corresponding dual events. A typical 
example is the event that the cluster of the origin does not touch the boundary dA, 
which is equivalent to the existence of a dual cluster in B(A*) containing a closed 
loop which surrounds the origin. 
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5.2. The Upper Bound. 



Our upper bound is stated in terms of the finite- volume approximation f^?^ to 



^wir' ^^^ equation (4.21). As in the last subsection, we will often explicitly indicate 
the /9-dependence of the relevant quantities. 

Theorem 5.2. Let d = 2, 0<(3<oo and q > I. Then there exists a constant 
Ci{(3,q) < cxD such that 

f^?r>(L,M)<Ci(/3,Q) (rf,ee,/3*((i^-l)ei))'. (5.9) 

Proof. Let A denote the box A(L,M), see equation (4.16). By its definition (4.21), 
the connectivity function t^°^«(-^, M) is the probability, in the measure /Uwir,/3,A, of 
the event i?o,L(^) = {0 ^ Lii} H {0 ^ dA}. Equivalently, Roj^{A) can be defined 
as the intersection of the event {0 ^^ Lii} and the event that there is a closed loop 
7* of occupied dual bonds surrounding the points and Lei. Consider the points 

x^ = (-1/2, ±1/2) and y^ = (L + 1/2, ±1/2) 

in A*. Given the fact that the connection from to Lei niust occur without touching 
9A, it is clear that the dual loop 7* must consist of four pieces: the bond {x'*L, x'^), 
a path 7^ connecting the point x!j_ to the point y!^, the bond (yl^y^), and a path 
7* connecting the point y'^ to the point x'^. Moreover, the two paths 7^ : x^ — > y^ 
must occur in B{A*) \ {{x*_,x'^), {yl,y^)}. Let us denote by R*^, i?|j, RZ. and i?!j_ 
the four events described above, namely (dual) occupation of the bond (x!l,x!!j_), 
the bond {yl,y^) and some paths 7^ : Xj_ -^ y^. in B{A*) \ {{x'L^x'^), {yl,y^)}, 



+ ■ 



respectively. Then i?[J"^(A) = {0 ^ Lei} n Rl n R*j^n R*_ n R 

Consider now a configuration to G -R£ fl R'^ fl RZ. fl i?i!j_. It is an easily verified 
geometrical fact that uj G Ro\ (A) if and only if the dual cluster joining x!j_ to y^ 
and the dual cluster joining xl to y'^ are connected only via the two bonds {x'^^x'^) 
and (i/l, y'^), i.e. if and only if there is no dual connection between x'^ and x*^_ in the 
set B{A*) \ {(xl,x!j_), {yt,y^)}. Using finite energy in the form (2.24) to convert 
the two occupation events i?^ ^^^ Rr ^o the events that the bonds (a;l,x!j_) and 
{y-j y^) are vacant, and the duality relation (5.8) to transform the measure A*wir,/3,A 
into the free measure niree,p*A* , we therefore obtain 

f^°^/3(i^,M) < Ci{(3,q) ;Ufree,/3*A*({a;*_ ^ yl} n « ^ y*+} n {x*_ ^ <}) (5.10) 

with Ci{l3,q) < 00 if < P < 00. 

Next, we note that by Proposition 2.8, 

l^bee,p*A*{{X- ^ y*-} n {xX ^ y+} n {X*_ 4A X+}) 

< /Ufree,/3*A* (2;+ ^ 2/+) /^free,/3*A* (^^l ^ y*_) ■ (S-H) 

Using the monotonicity (2.18), we obtain 

/Ufree,/3*A* (a^± ^ ?/±) < /Ufree,/3* (a^± ^ y±) = Tfree,/?* ((-^ ± l)ei) , (5-12) 

which, combined with (5.10) and (5.11), proves the theorem. D 
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5.3. The Lower Bound. 

Our lower bound is given in terms of the approximation fj^^g to Tf[^'^, see equation 
(4.23). 

Theorem 5.3. Let d = 2, 0<(3<oo and q > I. Then there exists a constant 
C2(/3, q) > such that for all positive integers M and L 

rXp(L^i)>C2{(3.q)''{fte,f3'{L-l.M)f. (5.13) 

Proof. In order to prove the theorem, we introduce several sets in both Z^ and in 
its dual (Z*)^. Consider the dual boxes 

a; = {x* G (Z*)2 I i < X* < L - i, l<x*2<2M + ^}, 

A*_ = {x* e (Z*)2 I 1 < a:* < L - i, -(2M + i) < x* < -i}, 

dual points 

x*^ = iei ± (M + i)e2, x*± = {L- i)ei ± (M + i)e2, 

dual bonds 

^±,1 = {x*±,x*^-ei), 6±_^ = (y±, y± + ei), 

and dual vertical lines 7^* and 7* joining the point xZ. — ei to the point x!j_ — ei, 
and the point y!i + ei to the point y^ + ei, respectively. In addition to these sets 
in (Z*)^ and 182, consider the points 

x± — ±Me2 and y± = Lei ± Me2, 

in Z^ and the vertical lines 7; and 7^. in 182 that join the point X- to the point x+, 
and the point y_ to the point y+, respectively. 

Consider now the events i?j_ that the points x'^ and y^ are connected by a path 
of dual occupied bonds in Aj_, and the events i?* that all bonds in 7* U 6!j_ q, U 6!l ^ 
{a = l,r) are occupied. The event i?!j_ fl i?l fl R^ fl i?* then clearly implies the 
existence of an occupied dual path surrounding the points and Lei, so that the 
event {0 ^ Lei} fl i?!j_ n i?!L fl i?f fl R* is contained in the desired event Rq^^ = 
{0 ^^ Lei} n {|C(0)| < 00}. As a consequence 

rfi[;,;3(^ei) = ;Uwir,/3(i?o,L) > /"wir,/3 (^{0 ^ Lci} H R^ H R*_ H Rl n i?;") . (5.14) 

Our goal is to modify the event in the argument of (5.14) so that (1) the event 
{0 *-^ Lei} is guaranteed to occur, and (2) the two dual paths across AZ^ and 
K*_ carry decoupling events that allow us to apply Proposition 2.6 to factor their 
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probabilities. We begin by degrading our estimate (5.14) by constructing vertical 
lines of occupied (direct) bonds: 

T^l^f^iLei) > ^i„ir,f3 Mc^7iU7. = 1} n {0 ^ Lci} H R^ H R*_ H R^ n R^Y 

Here, as usual, ujb denotes the configuration uj restricted to B. Using finite energy 
in the form (2.24) to fiip the 4M + 6 bonds in 7f U b*_^i U b*_i U 7; U b^^^ U b*_^^, we 
then obtain 

rX^iLei) > C{(3, Q)^^+Vwir,/3 (w^^u^^ = 1} n {0 ^ Lei} nR^HR: 

with a suitable constant C(/3, q) > 0. 

Consider now the events (i?^)^" that x^. and y^. are connected by dual cluster 
C*{xjj.) C B{Aj^), i.e. by clusters that lie entirely within i?(A^), and hence are 
surrounded by decoupling circuits of occupied (direct) bonds in (i?"'"(A^)) . Clearly 
Rl D (i?i)fi" and thus 

rX^iLei) > C{P, ?)^^+Vwir,/. ({u;7.U7. = 1} n {0 ^ Lei} n {R*+f'' n {R*_f-^ . 

(5.15) 
We now claim that 

{u^.uj. = 1} n {0 ^ Lei} n (i?;)fi" n {R*_f^ = {u^,u^^ = 1} n {R^f^ n {R*_f^ . 

(5.16) 
In order to see this, let a; be a configuration in {R"^)^"^. As noted above, the 
condition C*(x^) C B(A*) implies the existence of a closed path of occupied bonds 
in (i?+(A;!j_)) surrounding x^ and y^. Given u> E (-R!J_)*^", let 7 be the innermost 
such path. Since 7 surrounds x"^ and y^, but lies within (i?+(A;!j_)) , it must visit 
the points x^ and y_|_; thus it provides a connection between x^ and y+ by a path of 
occupied bonds. Observing that the vertical paths 7^ and 7^ connect the point to 
the point x+ and the point Lei to the point y+, we see that there is automatically 
an occupied path from to Lei, which completes the proof of (5.16). 

Using once more finite energy, the relations (5.15) and (5.16) together with the 
duality relation (5.8) now imply that 

> C(/?,Q)S^+Vwir,/3((i?;)^" n {R*_f^) 

= Ci(3, qf^+^^liree,P' [R^ ^ ^-") , (5-17) 

where R^ are the events dual to (i?^)^". 
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Finally, we use the fact that R^ and R^ are events of the form considered in 
Proposition 2.7, so that 

//free,/?* {Rt" n i?fi") > //free,/3* {R^^ /"free,/?* (i?^") = (r^^e (^, ^))'- (5-18) 

This completes the proof of Theorem 5.3. D 

Notice that, in contrast to the proof of Theorem 5.2, the above proof does not 
invoke monotonicity properties which depend on boundary conditions. Thus it can 
be used equally well to give a lower bound on t^^^, namely: 

Corollary. Let d = 2, 0<(3<oo and q > I. Then there exists a constant 
C2(/3, q) > such that for all positive integers M and L 

T|e\/3(^ei) >C2(/3,Q)^(f^°>(L-l,M))'. (5.19) 

5.4. The Dichotomy and Percolation Probabilities. 

In this subsection we prove Theorem 5.1. We start with a proposition which is 
essentially a corollary to the upper and lower bounds of Theorems 5.2 and 5.3. 

Proposition 5.4. Let d = 2,0<(3<oo and q > 1. Then 

CliP) = Uiree{n ^f P'^iP*) = (5.20) 

and 

^wfr(/3)=6ree(/3) ^/ P^ iP) = 0- (5-21) 

Proof. Introducing the notation $,^"^{M]/3) and ^|^gg(M;/3) to indicate the /3-de- 
pendence of the correlation lengths ^^°^(M) and ^^ee(-^) corresponding to f^?'^ 
and fjljee , the upper and lower bounds of Theorem 5.2 and 5.3 imply that 

|^°f(M;/3)<i^free(r) (5.22) 

and 

lilTeiM;n<^^M (5.23) 

Taking the limit M — > oo, and observing that the left hand side of (5.22) goes to 
^wfrl/^) t>y Lemma 4.8 and the corollary at the end of Section 4, while the left hand 
side of (5.23) goes to Ch^eif^*) t>y the same corollary and Lemma 4.7, we conclude 
that 

lutein < eir(/3) < |6ree(/3*). (5.24) 

Since ^f^eif^*) = ^free(/3*) if P^^'^if^*) = 0, this implies the first part of the proposi- 
tion. If PZ''{/3) = 0, then CUP) = ^wir(/9). In addition, by the results of [ACCN], 
A^wir,/? = A*free,/3 whenever P^"'(/?) = and hence CbeeiP) = ^wir(/?)- This implies 
the second part of the proposition. D 

In order to complete the proof of Theorem 5.1, we use the fact, proven in Section 
2.4, that the free measure /Ufree is ergodic under any nontrivial subgroup of the 
translation group (Theorem 2.10). Since, in addition, ^Ufree is an FKG measure 
which is invariant under horizontal and vertical translations and axis reflections, a 
bond percolation analogue of the theorem of [GKR] applies, leading to 
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Theorem 5.5. Let d = 2,0</3<oo and q > 1, and assume that P^^^{/3) > 0. 
Then, with probability one with respect to the free measure //free,/?, any finite set of 
sites in 1? is surrounded by a circuit of occupied bonds. 

Corollaries. Let d — 2,0<l3<oo and q >1. Then 
1) 

P'^'''mPZ''in = 0- (5.25) 

2) P^^{(3) = for all (3 > (3o ■ In particular, P^^'^{(3) is left continuous at (3o- 

3) If P^'{I3) = or P^''^(/3) > 0, then /Ufree,/3(-) = /Uwir,/3(-). o,nd in particular. 

Proof. As noted above, the theorem follows from (a bond percolation analogue of) 
[GKR] and Theorem 2.10. Corollary 1 then follows immediately, and Corollary 2 
also follows easily - see equation (1.14) and the paragraph preceding it. The first 
part of Corollary 3, namely that P^^{(3) = implies Hiree,i3i-) = A*wir,/3(')7 is a result 
of [ACCN]. That equality of the measures is also implied by P^^'^{(3) > follows 
from (5.25), [ACCN] and the self-duality of the model. D 

Remarks. 

1) Theorem 5.1 (the dichotomy) now follows immediately from Proposition 5.4 and 
equation (5.25). 

2) We expect that left continuity of P^^{(3) at the transition point holds in all 
dimensions provided q > 1. However, we do not expect P^^'^{(3) to be right 
continuous at the transition point if the system has a first-order transition; in- 
deed, for q sufficiently large, this can be established using Pirogov-Sinai theory, 
as used e.g. in [LMMRS]. This is to be contrasted with the behavior of P^^{(3). 
By standard percolation arguments [R], namely expressing P^^{(3) as the de- 
creasing limit of the finite-volume quantities (3.4) (which are continuous and 
non-decreasing in /3), P^"^(/9) is right continuous for all (3 and all g > 1 in di- 
mension d > 1. However, in dimension d > 2, convergent expansions have been 
used to show P^'^{(3) is not left continuous at the transition point provided q is 
sufficiently large [KoS] (see also [LMR] and [LMMRS]). 

3) Corollary 3 implies that in two dimensions the Gibbs state is unique at all /3 
except those for which P^'^((3) = while P^^{(3) > 0. Presumably, this never 
occurs for systems with second-order transitions {q < 4md = 2), and occurs only 
at a single point - the transition point - for systems with first-order transitions 
{q > 4 in d = 2). Again, this can be proven via expansion methods in d > 2 for 
q sufficiently large. 

We conclude this section with a little result which is an easy consequence of 
Proposition 5.4. The result shows that continuity of the magnetization at (3o ensures 
criticality of the transition, i.e. divergence of the correlation length(s). 
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Proposition 5.6. Let d = 2 and q > 1. Then M{Po) = P^''{Po) = implies 

^wh-(/5o) = OO and hence also ^|ee(/5o) = 6ree(/3o) = ^wir(/3o) = oo. 

Proof. By the assumption P^^{(3o) = 0, the FKG ordering of states (2.23) and the 
definition of /9o, we have = P^'^{(3o) = P^'^iPo)^ ^^^ hence by the first branch 
(5.20) of Proposition 5.4: 

CliPo) = U^reem = '^Ciree{Po) • (5.26) 

On the other hand, again by the assumption P'^^{l3o) = 0, we have the second 
branch (5.21) of Proposition 5.4, namely: 

^lWo)=^ireem. (5.27) 

From (5.26) and (5.27), we conclude that either ^free(/3o) = or ^free(/3o) = oo. The 
first case is easily ruled out by considering e.g. rfi.ee(ei). That the other correlation 
lengths also diverge is an immediate consequence of remark following Theorem 
4.3. D 



Appendix: Reflection Positivity and the Transfer Matrix 

The concept of refiection positivity and its consequences are well-known tools in 
the context of field theory. For the convenience of the reader we give a brief review 
in this appendix. 

We consider a (finite or infinite) lattice A C Z'^ which is invariant under refiections 
at a plane E. Here S is either a lattice plane or a plane which lies halfway between 
two lattice planes. Denoting the refiection at E by r, we then decompose A as 
A = A+ U A_, where A+ are the points on one side of E, A_ = r(A+) are the points 
on the other side of E, and A_ fl A_|_ = E fl A (which is of course empty if E lies 
between two lattice planes). 

For a local observable A with support suppA C A_|_, one introduces the refiected 
observable r{A) as 

(r(A))(a) = A(r(a)) (A.l) 

where r{a)j: = ar(x)- Refiection positivity of the Potts model is the statement that 



{riA)A),^^>0. (A.2) 

The proof of (A.2) is standard; for the strategy, see e.g. [FILS], [Se]. 

The inequality (A.2) has several important consequences. Here we are mainly 
interested in the representation of truncated expectation values as matrix elements 
of a suitably defined transfer matrix T. In order to define the transfer matrix T 
in the setting considered here, we need, in addition to the refiection invariance of 
A, that A is invariant under translations perpendicular to E. We therefore assume 
that A is of the form 

A = Z X Ai 
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where Ai is a (finite or infinite) sublattice of 7/~^ . We then consider the algebra 
Aj^ of local observables with support in A+, where from now on 

A+ == {(x, f) G Z'^lx > 0, X G Ai} . 

while E = {(x, x) G Z'^ja; == 0, a^ G Ai}. Due to (A. 2), the equation 



{A,B):={r{A)B),^^ (A.3) 

defines a positive semi-definite scalar product over Aj^. Dividing out the corre- 
sponding null space N and completing the resulting space in the usual way, this 
leads to the definition of a Hilbert space Ti = A+/Af. 

Next, we introduce, for each local observable A G A.^-, the observable TA which 
is obtained from A by translation by one lattice unit in the positive direction per- 
pendicular to S. It is an easy consequence of the Cauchy-Schwarz inequality for the 
scalar product (A.3), see [Se] for details, that T obeys the inequalities 

< (A, TA) < (A, A) . (A.4) 

The operator T therefore defines a positive transfer matrix, which obeys the in- 
equalities 

0<T< 1 

as an operator on 7i. Observing that the vector O corresponding to the constant 
function 1 G A+ is an eigenvector of T with eigenvalue 1, we note that the norm of 
T is one. 

We finally consider the interpretation of truncated expectation values in the 
above Hilbert space representation. Since 



{riA)),^A = {riA)A),^^ = {An) (A.5a) 

while 

{T^A),^^ = (O, T^A) = (T"0, A) = (O, A) , (A.5b) 

one immediately obtains 



(r(A); T^A),^^ = (A, T^A) - (A, O) (O, A) . (A.6) 

Introducing the projection operator P± onto the Hilbert space orthogonal to O, 
equation (A.6) becomes 



{r{Ay,T-A),^^ = {A^,T-A^), (A.7) 

where 

A± = P±A = A - (O, A)0 . (A.8) 
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If the support of A is a subset the lattice plane E, r{A) = A, and equation (A. 7) 
reduces to 

{A;T-A),^^ = {A^,T-A^). (A.9) 



Equation (A.9) is an important technical tool in the proof of the existence of the 
correlation lenghts P -I and P -1. 

Remark:. In the context of Euclidean field theory, the direction perpendicular to E 
is often interpretated as the Euclidean time. The Hilbert space 7i = A+fN" is then 
nothing but the quantum mechanical Hilbert space of the considered model, and T 
is the generator of the Euclidean time translations, i.e. T = e"*^^, where e is the 
lattice spacing and H is the Hamilton operator of the theory. 

However, 7i and T have no such interpretation for the classical Potts model. This 
is due to the fact that here A is the lattice of a classical system, and not a lattice 
approximation to Euclidean space-time. 
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